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1.0  OBJECTIVE 


This  report  assesses  the  National  Aeronautics  and  Space  Administration  (NASA)  long-term, 
orbital-debris  engineering  model  developed  to  predict  the  debris  environment  in  low  earth  orbit 
(LEO).  Specifically,  this  report  is  an  assessment  of  the  NASA  orbital-debris  model  as  described  in 
References  1  and  2,  hereafter  known  as  the  NASA  89  and  NASA  90  models,  respectively.  This 
report  describes  the  model  development,  presents  a  model  sensitivity  analysis,  arid  establishes  the 
model  uncertainties.  The  goal  of  this  assessment  is  to  assist  the  Air  Force  (AF)  in  its  task  of 
developing  a  long-term  debris  model  that  can  (1)  operate  with  higher  accuracy  at  the  relevant 
altitudes  and  orbits  parameters;  (2)  benefit  from  new  AF  and  non-AF  debris  measurements;  and  (3) 
accommodate  current  and  future  space  scenarios. 


Section  2.0  presents  a  brief  background  and  overview  of  orbital  debris  and  why  this  environment 
is  of  concern  to  space  operators,  especially  the  Department  of  Defense  (DoD)  and  the  AF.  The 
NASA  89  model  and  the  revised  NASA  90  are  discussed  in  Section  3.0.  The  assessment  focuses 
on  the  revised  NASA  90  model.  The  NASA  89/90  model  development  is  presented  in  Section  4.0. 
Included  is  an  examination  of  the  measurement  data  bases,  model  development  rationale,  and  an 
assessment  of  the  model  limitations.  The  orbits  of  interest  (altitudes  and  inclinations)  for  current 
and  future  AF  and  DoD  operations  are  discussed  in  Section  5.0.  The  sensitivity  analysis  of  the 
NASA  90  model  is  found  in  Section  6.0.  The  model  sensitivities  are  analyzed  with  respect  to  the 
model  variables:  debris  diameter,  orbital  altitude  and  inclination,  time,  solar  activity,  and  debris 
growth.  Summary  comparisons  and  an  assessment  of  the  sensitivity  are  presented.  The  ability  of 
the  model  to  predict  the  environment,  and  its  associated  error  or  uncertainty,  is  covered  in 
Section  7.0.  This  uncertainty  analysis  for  the  NASA  90  model  is  based  upon  a  propagation  of 
error  analysis  that  uses  NASA's  estimated  uncertainties  to  predict  the  overall  model  uncertainty. 
Section  8.0  summarizes  the  assessment  of  the  NASA  90  model. 


1 


2.0  BACKGROUND 


The  man-made,  orbital-debris  environment  is  of  increasing  concern  to  users  of  space.  This 
environment  has  grown  since  man's  first  steps  into  space.  If  unchecked,  it  will  pose  extreme 
hazards  for  future  space  operations.  An  essential  first  step  in  assessing  potential  hazards  posed  by 
the  environment  is  to  characterize  the  actual  orbital-debris  environment  A  diagnostic  model  should 
be  based  on  complete  and  reliable  observational  data.  From  diagnostic  models,  predictive  models 
can  be  developed  that  may  be  based  on  either  physics  or  derived  empirically  from  the  measurement 
data.  The  NASA  90  model  is  an  empirically  derived,  predictive  model  of  the  debris  environment 
for  space  system  designers,  operators,  and  other  users.  The  emphasis  of  NASA  90  model  is  on 
the  prediction  of  the  long-term,  orbital-debris  environment  Before  the  NASA  90  model  can  be 
used  with  confidence,  it  must  be  assessed  against  the  actual  data  or  validated  models  to  determine 
its  limits  of  applicability  and  accuracy.  The  purpose  of  this  effort  is  to  conduct  such  an 
assessment 

Early  in  man's  first  operations  in  space,  the  principal  concern  with  possible  particle  impacts  was 
due  to  the  natural,  meteoroid  environment  Meteoroids  are  part  of  the  interplanetary  environment 
and  are  encountered  as  the  earth  sweeps  through  space  in  its  orbit  about  the  sun  (meteoroids  are  not 
in  earth  orbit).  At  any  one  instant,  it  is  estimated  that  the  meteoroid  mass  within  2(X)0  km  of  the 
earth's  surface  is  ~200  kg  with  most  of  this  mass  concentrated  in  0.1-mm-dia  meteoroids. 
Meteoroids  move  at  average  relative  speeds  of  20  km/s  (Ref.  1).  The  natural  environment  is  of 
concern,  but  differs  from  the  orbital-debris  environment.  The  man-made  orbital  debris  remains  in 
earth  orbit  and  originates  from  the  manner  of  space  operations  (intentionally  or  unintentionally 
leaving  debris,  boosters,  spent  payloads,  etc.,  in  orbit),  intentional  or  unintentional  explosions, 
and  collisions.  Estimates  of  orbital  debris  now  in  earth  orbit  within  2000-km  orbital  altitude  ranges 
anywhere  from  1.6  to  3.0  Mkg.  Orbital  debris  moves  at  relative  speeds  on  the  order  of  10  km/s. 
Using  the  value  3.0  Mkg  for  orbital-debris  mass  in  orbit,  it  is  further  estimated  that  -300  kg  of 
this  mass  consists  of  particles  <1  mm,  1000  kg  consists  of  particles  <1  cm,  with  the  remainder 
mostly  concentrated  in  approximately  3000  spent  rockets  and  inactive  payloads  (Ref.  1).  Unlike 
the  meteoroid  environment,  much  of  the  orbital  debris  is  distributed  in  regions  of  current  satellite 
operations.  The  orbital-debris  environment  has  become  at  least  as  hazardous  as  the  meteoroid 
environment  for  objects  in  LEO. 


3.0  NASA  ORBITAL-DEBRIS  MODEL 


This  section  describes  the  NASA  89  and  90  orbital-debris  models.  Included  in  these  models  is  the 
prediction  of  (1)  the  orbital-debris  flux  (impacts  per  year  per  unit  area  (square  meters))  as  a 
function  of  the  debris  diameter,  year  of  interest,  solar  activity,  spacecraft  altitude  and  orbital 
inclination,  and  projected  growth  of  the  orbital  debris;  (2)  the  orbital-debris  collision- velocity 
impact  and  direction  distribution  as  a  function  of  orbital  inclination;  and  (3)  the  orbital-debris  magg 
as  a  function  of  diameter  for  assumed  spherical  shapes. 

The  model  was  developed  with  the  goal  of  producing  ah  engineering  model  for  use  by  the 
engineering  and  design  commumty  to  apply  in  their  regions  of  interest.  The  model  represents  a 
curve  fit  to  the  measurement  data  (data  sources  and  model  development  will  be  discussed  in 
Section  4.0).  The  model  is  considered  to  be  2q)plicable  only  to  objects  in  LEO  (Refs.  1  and  2). 
This  represents  a  nominal  altitude  of  h  ^  2000  km;  however,  Kessler*  felt  it  was  most  valid  for 
altitudes  of  h  ^  1000  km.  Subsection  3.4  examines  the  behavior  of  the  model  in  the  extreme  limits 
of  the  variables  on  which  the  model  is  based. 


3.1  NASA  89  ORBITAL-DEBRIS  R.UX  MODEL 


The  NASA  89  orbital-debris  flux  engineering  model  (Ref.  1)  predicts  the  cumulative  flux  of  orbital 
debris,  for  a  diameter  (d)  and  larger,  for  space  systems  in  LEO  at  altitude  (h),  inclination  (i),  and 
year  (t)  using  the  following  expression: 


F(d,  h,  i,  t,  S)  =  k  (h,  S)  ^  (i)  [Fi(d)gi(t)  +  F2(d)g2(t)]  (1) 

where 


F 

d 

h 

i 

t 

S 


=  Time-averaged,  cumulative,  orbital-debris  flux,  impacts  per  year,  and  per  unit 
area  (impact/yr  •  m  2)  against  a  single-sided  surface  for  diameters  (^d) 

=  Orbital-debris  diameter  (cm) 

=  Altitude  (km) 

=  Orbital  inclination  (deg) 

=  Time  (yr) 

=  Thirteen-month,  smoothed,  solar  radio  flux  at  10.7-cm  wavelength  expressed 
in  10  4  Jansky  (Jy)  for  the  previous  year  (t  -  1),  typically  70  ^  S  ^  150  where 
S  may  exceed  value  of  200  at  solar  maximum 


* 


Private  meetings/communication  with  NASA/JSC,  D,J.  Kessler,  1990-1991. 
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k 


‘J>(h,S) 


^i(h.S) 

'F(i) 


Fi(d) 

gi(t) 

P 

F2(d) 

g2W 


Surface  orientation  factor  where  k  is  the  ratio  of  the  flux  on  a  single-sided 
surface  at  a  specified  orientation  to  the  flux  on  a  randomly  tumbling  single¬ 
sided  surface,  k  =  1  for  a  randomly  tumbling  surface  and  theoretically, 

0  <  k  <  4,  for  a  fixed  surface  orientation  (Ref.  1) 

Oi(h,S)/[<I>i(h,  S)  +  1],  functional  form  accounting  for  effects  of  solar 
activity  (atmospheric  drag  effects  on  debris)  as  a  function  of  altitude  and  the 

solar  radio  flux 
10  (h/200  -  s/140  -  1.5) 

Debris  concentration  correction  factor  (Table  1)  where  4^(i)  is  the  ratio  of 

the  flux  at  inclination  i  to  the  flux  at  the  population's  average  inclination  of 

approximately  60  deg  or  Flj  /  F||  _  go  deg 

1.05  (10  -5)  d  -2-5,  size  distribution  for  small  particles 

(1  +  2p)  (t-  1985)^  small-particle  growth 

Assumed  annual  growth  rate  of  mass  in  orbit  (5%  or  p  =  0.05) 

7  (10  10)  (d  -t-  7(X))  size  distribution  for  large  particles 
(1  -I-  p)  (H985)^  large-particle  growth 


Note  that  the  distinction  between  small  and  large  orbital-debris  particles  is  in  the  range  of  1  to 
10  cm.  No  distinct  boundary  is  specified  in  the  NASA  reports;  however,  diameters  <1  cm  are 
considered  small  by  the  technical  community  while  diameters  >10  cm  are  considered  large. 

The  average  number  of  impacts,  N,  on  a  surface  area.  A,  exposed  to  the  orbital-debris 
environment,  F,  over  a  period  of  time,  t,  to  tf,  is  calculated  using  the  following  integral 
relationship: 


N  =  FAdt  (2) 

Note  that  Equation  2  determines  the  average  number  of  impacts  on  a  single-sided  surface  whose 
orientation  is  given  by  the  value  of  k.  To  find  the  total  number  of  impacts  on  a  spacecraft,  one 
must  determine  k  for  each  respective  surface  and  then  sum  up  the  impacts  evaluated  for  each 
surface,  e.g.. 
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for  M  surfaces.  The  only  value  changing  in  the  orbital  flux  expression,  (Eq.  1)  is  the  value  of  k, 
therefore 


M  ( 

^TOTAL  =  X  ~ |j 


Table  1.  Orbital-debris  concentration  correction  factor,  'F  (i). 


Inclination  'F(i)  Inclination  'F(i)  Inclination  'F(i) 

(deg)  (deg)  (deg) 


25 

0.900 

58 

26 

0.905 

59 

27 

0.910 

60 

28 

0.912 

61 

28.5 

0.9135 

62 

29 

0.915 

63 

30 

0.920 

64 

31 

0.922 

65 

32 

0.927 

66 

33 

0.930 

67 

34 

0.935 

68 

35 

0.940 

69 

36 

0.945 

70 

37 

0.950 

71 

38 

0.952 

72 

39 

0.957 

73 

40 

0.960 

74 

41 

0.967 

75 

42 

0.972 

76 

43 

0.977 

77 

44 

0.982 

78 

45 

0.990 

79 

46 

0.995 

80 

47 

1.000 

81 

48 

1.005 

82 

49 

1.010 

83 

50 

1.020 

84 

51 

1.025 

85 

52 

1.030 

86 

53 

1.040 

87 

54 

1.045 

88 

55 

1.050 

89 

56 

1.060 

90 

57 

1.065 

91 

1.075 

92 

1.400 

1.080 

93 

1.440 

1.090 

94 

1.500 

1.100 

95 

1.550 

1.115 

96 

1.640 

1.130 

97 

1.700 

1.140 

98 

1.750 

1.160 

99 

1.770 

1.180 

100 

1.780 

1.200 

101 

1.770 

L220 

102 

1.750 

1.240 

103 

1.720 

1.260 

104 

1.690 

1.290 

105 

1.660 

1.310 

106 

1.610 

1.340 

107 

1.560 

1.380 

108 

1.510 

1.410 

109 

1.460 

1.500 

no 

1.410 

1.630 

111 

1.380 

1.680 

112 

1.350 

1.700 

113 

1.320 

1.710 

114 

1.300 

1.700 

115 

1.280 

1.680 

116 

1.260 

1.610 

117 

1.240 

1.530 

118 

1.220 

1.490 

119 

1.200 

1.450 

120 

1.180 

1,410 

121 

1.165 

1.390 

122 

1.155 

1.380 

123 

1.140 

1.370 

124 

1.125 

1.380 

125 

1.110 
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Finally,  the  probability  of  exactly  n  impacts  occurring  on  a  single-sided  surface  is  found  using 
Poisson  statistics  where 


3.2  NASA  90  ORBITAL-DEBRIS  FLTJX  MODFT. 

The  NASA  90  orbital-debris  flux  engineering  model  (Ref.  2)  is  an  update  to  Equation  1  and  is 
based  upon  more  recent  measurements  of  the  orbital-debris  environment  These  updates  will  be 
discussed  later  in  this  report.  The  revised  model  is  as  follows: 

F(d.  h,  i,  t,  S)  =  H(d)  a)(h,S)  Td)  [Fi(d)  gi(t)  F2(d)  g2(t)]  (6) 


where 


F(d,h,i,t,S)  = 

Time-averaged,  cumulative,  orbital-debris  flux  on  a  randomly 
tumbling  single-sided  surface  for  diameters  d  and  greater.  Note  that 
the  parameter  k  is  no  longer  associated  with  F,  but  is  now  only 
associated  with  the  calculation  of  the  average  number  of  impacts,  N, 
on  a  surface  of  A,  defined  in  Equation  7. 

H(d) 

f  10«^P(-(  log  lod- 0.78)  V(  0.637)2) 

Fi(d) 

1.22  (10  -5)  d  -2-5 

F2(d) 

8.1  (10  10)  (d -H  700) -6 

gi(t) 

(1  +  q)  (1  - 1988)  where  q  =  0.02  until  2010  and  q  =  0.04  after  2010 

g2(t) 

1  -1-  p  (t  -  1988)  where  p  =  0.05 

All  other  parameters  and  functions  remain  as  previously  defined. 


N  = 


M  \ 

IkjFAjdtj 


(7) 


The  probability  of  exactly  n  impacts  remains  as  defined  in  Equation  5. 
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A  comparison  of  the  NASA  89  model,  (Eq.  1)  to  the  NASA  90  model  (Eq.  6),  is  shown  in  Figure 

1  for  h  =  500  km,  t  =  1995,  i  =  30  deg,  S  =  90,  and  k  =  1.  The  principal  difference  is  due  to  the 
function  H(d)  in  the  region  1  cm  <  d  <  100  cm.  In  Subsection  4.3.3,  H(d)  will  be  shown  to 
increase  the  NASA  89  prediction  for  these  large  particles  by  up  to  a  factor  of  3.  Other  differences 
in  the  model  flux  predictions  between  the  NASA  89  and  90  models  are  due  to  changes  in  the 
growth  of  both  large  and  small  debris.  The  NASA  90  model  assumed  growth  and  growth  rates  are 
less  severe.  For  the  small  debris,  the  growth  rate  has  been  reduced  from  10  percent  (2  x  p)  to 

2  percent.  The  reduction  in  growth  rate  is  based  on  studies  that  use  tlie  NASA  EVOLVE 
(Evolutionary  Model)  code  (Ref.  2).  EVOLVE  shows  fewer  events  and  a  lower  source  term  from 
collisions  of  larger  objects.  This  reduction  in  the  predicted  environment  for  the  small  debris  can  be 
seen  in  Figure  1.  For  the  large  debris,  the  functional  form  of  the  large-particle  growth,  g2(t),  has 
been  changed  from  compounded  to  linear.  The  reason  for  this  will  be  discussed  in  Section  4.0. 
The  assumed  growth  rate,  p,  remains  at  5  percent.  Comparisons  of  the  NASA  90  model 
predictions  versus  the  measured  environment  are  shown  in  Figures  2  and  3. 


Diameter  (cm) 

Figure  1.  NASA  89  model  versus  NASA  90  model,  flux  versus  diameter, 
h  =  500  km,  i  =  30  deg,  t  =  1995,  S  =  90,  (k  =  1). 
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Figure  2.  NASA  90  orbital-debris  flux  versus  diameter,  d  (cm),  compared  to  data  from 
Solar  Max,  Goldstone,  Aricebo,  and  GEODSS  (h  =  500  km,  i  =  28.5  deg, 
t  =  1988,  S  =  140). 


Figure  3.  NASA  90  orbital-debris  flux  versus  altitude,  h  (km),  compared  to 
USSPACECOM  data  (d  >  10  cm,  i  =  60  deg,  t  =  1990,  S  =  200). 
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3.3  DEBRIS  IMPACT  AND  DIRECTION  DrSTRTBTmONS.  AND  PARTICLE  DENSTTY 

The  orbital-debris,  collision- velocity  impact  distribution  and  impact  direction  remains  unchanged 
between  the  NASA  89  and  90  models.  The  collision-velocity  distribution  model  predicts  the 
number  of  impacts  with  collision  velocities  between  V  and  V  +  dV  relative  to  a  spacecraft  as  a 
function  of  the  spacecraft's  orbital  inclination,  i,  and  the  collision  velocity,  V  (km/s),  where 


A  is  a  constant  and  B,  C,  D,  E,  F,  G,  H,  and  Vq  are  functions  of  the  spacecraft  orbital  inclination. 
The  parameters  A  through  H  result  from  the  empirical  curve  fit  while  V©  is  related  to  the  debris 
average  orbital  velocity  as  seen  from  the  data  for  these  inclinations.  (Note  that  NASA  has 
apparently  introduced  the  appropriate  units,  i.e.,  impacts  •  s  2/km  2,  into  the  coefficients  to  give 
f  (V)  the  proper  units,  impacts.) 


f(V)  =  Number  of  impacts  between  V  and  V  -t-  dV,  where  f(V)  >  0  (if  a  value  of  f(V)  <  0 
results  from  Equation  8,  then  set  f(V)  =  0) 


V  _  /  7.25  +  0.015  (i  -  30) 

°  *■  1  7.7 

A  =  2.5 

-  (  0J-0.1(i-60) 


i  <  60  deg 
i  >  60  deg 


i  <  60  deg, 

60  <  i  <  80  deg 
i  >  80  deg 


^  f  0.0125  i<  100  deg 

^  =  1  0.0125  +  0.00125  (i  - 100)  i  >  100  deg 

D  =  1.3  -  0.01  (i  -  30) 

E  =  0.55  -I-  0.005  (i  -  30) 


F 


0.3  0.0008{  i  -  50)  2 
0.3  -  0.01  (i  -  50) 

0.0 


i  <  50  deg 
50  <  i  <  80  deg 
i  >  80  deg 
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G  = 


18.7 

18.7  +  0.0289  (i-60)3 
250 


i  <  60  deg 
60  <  i  <  80  deg 
i  >  80  deg 


H  = 


1.0  -  7.57  (10  -5)  (i  -  60)  2 


f(V)  represents  an  averaged,  collision-velocity  impact  distribution  over  all  altitudes  (LEO). 
Figure  4  shows  f(V)  for  an  orbital  inclination  of  30  deg. 


Figure  4.  Orbital-debris,  collision- velocity  impact  distribution, 
f(V),  versus  velocity  (V)  for  i  =  30  deg. 


The  normalized  velocity  distribution  is  given  by  (Refs.  1  and  2) 


ftV)=f{V)/  f(v)dV  (9) 

J  o 

Only  values  of  f(V)  >  0  have  any  physical  meaning,  therefore,  the  region  of  f(V)  that  is  used  to 
specify  the  collision  velocity  impact  distribution  is  given  from  f(V)  =  0  at  V  =  0,  defined  as  Vmin. 
to  f(V)  =  0  at  V  =  Vmax-  Equation  10  shows  that  the  positive  area  of  the  curve,  f(V),  is  used  to 
normalize  the  velocity  distribution.  This  is  given  as. 


f'{V)=f{V)/  f{V)dV 


i 


max 


(10) 
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where  Vmin  and  Vmax  represent  the  minimum  and  maximum  velocities  at  f(V)  equals  zero.  For 
each  orbital  inclination,  i,  Vmax  has  some  value  that  is  a  function  of  i.  The  maximum  Vmax , 

V'max.  for  all  inclinations  is  15.4  km/s. 

Once  the  collision-velocity  impact  distribution  and  velocities  are  known  for  a  given  orbital 
inclination,  then  the  direction  of  impact  as  well  as  the  frequency  of  impact  for  a  given  direction  may 
be  estimated.  This  is  done  by  first  assuming  the  direction  of  impact  is  specified  by  the  angle  a 
between  the  spacecraft  orbit  and  relative  impact  vector  as  shown  in  Figure  5.  The  spacecraft  and 
debris  velocity  vectors  are  assumed  to  be  circular  orbit  velocity  vectors  in  the  local  horizontal  plane 
B.  The  relative  velocity  is  the  resultant  vector  between  the  spacecraft  and  debris  inpitial  vectors. 


z 


Figure  5.  Orbital-debris/spacecraft  reference  frame. 
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From  the  data,  an  average  orbital  velocity  for  objects  in  LEO  of  7.7  km/s  was  used  to  determine  the 
collision  direction.  (Note  that  7.7  km/s  corresponds  to  a  circular  orbit  of  345  km,  but  the  higher 
velocity  more  accurately  reflects  the  existence  of  debris  in  elliptical  orbits.)  Using  the  geometry 
presented  in  Figure  5,  it  is  clear  that  the  maximum  impact  velocity  is  when  a  =  0  (head-on 
collision)  and  diminishes  to  0  when  a  >  90  deg.  A  cosine  relationship  provides  the  direction  angle 
of  impact  (a)  given  an  impact  velocity,  V,  and  the  maximum  possible  impact  velocity,  V'max  = 

15.4  km/s  (i.e.,  2  x  7.7  km/s). 

cos  (±a)  =  V  /  V'max  (11) 


The  range  of  values  for  V  is  bounded  by  Vmitt  and  Vmax  for  a  given  inclination. 


The  assumption  of  circular  orbits  implies  no  out-of-plane  (pitch  angle)  impacts  are  considered  in 
this  model.  This  is  an  approximation  of  the  actual  environment  since  debris  does  exist  in  elliptical 
orbits  and  will  have  some  small  pitch  angle  relative  to  the  local  horizontal  plane  containing  the 
spacecraft  velocity  vector.  This  approximation  is  valid  if  it  is  assumed  that  the  mostly  near  circular 
eccentricity  distributions  observed  for  large  trackable  population  holds  for  all  debris  sizes. 


Finally,  the  average  mass  density  for  orbital  debris  is  given  as  follows: 


P  = 

P  = 


2.8d^'^'^  g/cm^,  d>l  cm 
2.8  g/cm^,  d<l  cm 


}  (NASA  89) 


2.8d^'^^  g/cm^,  d>0.5  cm 
4.7  g/cm^,  d<0.5  cm 


}  (NASA  90) 


(12) 


For  purposes  of  the  model,  the  orbital  debris  particles  are  assumed  to  be  spherical,  thus  the  mass  is 
given  as  follows: 


3.4  OBSERVATIONS  rLIMITTNO  VATTTF..^;^ 


The  NASA  orbital-debris  model  is  an  engineering  tool  that  is  based  on  a  curve  fit  to  data.  The 
model  does  not  necessarily  follow  from  physical  laws,  therefore,  should  not  be  expected  to  be 
representative  in  regions  where  there  are  no  data  to  support  the  model.  One  assessment  of  the 
model's  performance  is  to  observe  its  behavior  when  its  variables  are  taken  to  the  limit.  Physical 
intuition  is  used  to  assess  the  model's  validity  at  those  limits;  this  will  establish  the  qualitative 
rairge  of  validity  and  will  point  out  any  regions  where  its  formulation  may  be  invalid.  This 
examination  is  not  meant  to  replace  other  analysis  necessary  to  assess  the  model.  Rewriting  the 
NASA  90  model  in  terms  of  its  variables  results  in  the  following: 


F(d,h,i,t,S) 


[(  -  log  10  d  -  0.78 )  2/(  0.637  )  2] 

=  V  10® 


X 


M 


(  h/200-S/140-1.5  ) 


X  VF(  i) 


V  10  {  h/200-S/140-1.5  )  ^2 

x([  1.22(  10*^)  d*2-^]  ( 

-K  8.1)  {  10^°)(  d  +  700)-^[  1-Hp(  t-  1988)] 


(14) 


The  limits  of  interest  can  be  observed  by  examining  the  behavior  of  F(d,  h,  i,  t,  S)  as  diameter  (d), 
altitude  (h),  and  growth  (p  and  q)  approach  both  zero  and  infinity  and  time  (t)  approaches  infinity. 
Note,  there  is  no  physical  significance  in  the  limiting  behavior  of 'F(i)  or  F  with  respect  to  orbital 
inclination,  i.  The  results  of  these  observations  are  as  follows: 


Limiting  Value 

As  d 

0, 

F 

— >  oo 

Observations 

Not  unreasonable 

As  d 

F 

0 

Reasonable 

As  h 

0, 

F 

Min  Constant  ^  0 

Not  unreasonable 

As  h  — > 

F 

Max  Constant 

Unreasonable 

As  d 

F 

— >  oo 

Reasonable 

As  p,q  — > 

0, 

F 

->  constant 

Reasonable 

As  p,q  — > 

F 

— >  oo 

Reasonable 
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As  diameter  decreases,  one  might  expect  the  number  of  particles  or  flux  to  increase  (more  smaller 
particles  than  larger  particles)  down  to  some  limiting  diameter;  therefore,  this  behavior  in  F  is  not 
unreasonable.  Likewise,  as  diameter  increases,  one  would  expect  the  flux  to  go  to  zero  simply 
based  on  physical  reality;  therefore,  on  this  basis  the  behavior  of  F  is  reasonable.  As  altitude 
decreases,  one  would  expect  (due  to  atmospheric  effects)  that  the  flux  environment  will  decrease, 
which  it  does.  In  the  limit,  h  =  0,  i.e.,  the  earth's  surface,  the  behavior  of  F  is  to  decrease  to  a 
minimum  value,  as  altitude  decreases,  which  seems  to  comply  with  physical  intuition  and  is  not 
unreasonable.  One  would  expect  that  F  should  decrease  as  altitude  increases  and  approach  zero  in 
the  limit  on  the  basis  of  physical  reality  (because  this  model  represents  earth  orbital  debris).  The 
behavior  of  F  as  altitude  increases  gives  evidence  that  the  model  is  valid  only  to  a  limiting  altitude. 
Its  use  should  be  restricted  to  LEO  as  specified  (Refs.  1  and  2).  The  behavior  of  F  as  time 
increases  is  reasonable  assuming  an  ideal  external  space  environment  (no  perturbations)  and  no 
change  in  growth  rates.  The  behavior  of  F  with  respect  to  growth  rates,  p  and  q,  is  reasonable 
especially  for  unconstrained  growth.  One  also  should  note  that  <I>(h,S)  behaves  as  expected  for 
increasing  or  decreasing  values  of  h  and  S  (as  will  be  shown  graphically  in  the  next  section). 
Qualitatively,  the  model  behaves  reasonably  in  these  limits  except  for  increasing  altitude. 
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4.0  MODEL  DEVELOPMENT/RATIONALE 


This  section  presents  the  rationale  and  basis  for  the  development  of  the  NASA  orbital-debris 
model.  As  mentioned  previously,  the  NASA  orbital-debris  model  is  an  engineering  model  and  is 
based  on  a  curve  fit  to  the  measurement  data  for  objects  in  LEO.  It  is  felt  by  its  author*  to  be  most 
accurate  for  h  ^  1000  km.  The  data  base  used  in  the  development  of  this  model  will  be  discussed 
as  weU  as  the  assumptions  and  development  of  the  functional  forms  resulting  from  the  curve  fit 
The  estimated  model  uncertainties  and  qualitative  limitations  are  also  described. 

4.1  DATA  BASF, 


The  NASA  orbital-debris  engineering  model  (Eq.  6)  and  the  collision-velocity  impact  distribution 
(Eq.  8)  represent  curve  fits  to  optical  and  radar  ground  observations  of  orbital  debris  and  data  from 
returned  spacecraft  surfaces.  Data  used  in  the  model  development  and  subsequent  revision  are  as 
follows: 

•  The  U.S.  Space  Command  (USSPACECOM)  radar  tracking  data  (Ref.  3),  both 
catalogued  and  uncataloged  Orbital  Element  Sets,  were  used  to  establish  altitude  and 
solar  activity  effects,  size  distribution  particles),  historical  growth  rates, 
inclination,  velocity,  and  directional  effects.  The  data  used  represent  "snapshots"  of  the 
LEO  environment  taken  once  a  year  from  1976  through  1988,  normally  in  February. 
These  data  provide  administrative  information,  period,  inclination,  apogee,  perigee,  and 
radar  cross  section  of  the  tracked  objects.  Sample  data  from  May  1990  are  shown  in 
Figure  6.  The  data  from  the  orbital  element  sets  have  been  converted  to  flux  units 
(impacts/yr  •  m  2).  The  USSPACECOM  currently  tracks  approximately  6500  to  7000 
objects  (Ref.  1). 

•  Optical  measurements  include  data  from  the  Massachusetts  Institute  of 
Technology/Lincoln  Laboratory  Experimental  Test  Site  (MIT/LL  ETS)  in  Socorro, 

New  Mexico  (Ref.  4)  and  the  Ground-based  Electro-Optical  Deep  Space  Surveillance 
(GEODSS)  network  (Ref.  5).  The  GEODSS  sites  are  located  near  Socorro,  as  well  as 
Maui,  Hawaii,  and  the  island  of  Diego  Garcia  in  the  Indian  Ocean.  The  GEODSS  data 
were  used  to  substantiate  altitude  effects  and  size  distribution.  The  MIT/LL  ETS  data 
consisted  of  10  hr  of  viewing  in  1984  and  10  hr  in  1985.  The  GEODSS  data  represent 
80.9  hr  of  viewing  (Ref.  5).  Sample  GEODSS  data  are  shown  in  Figure  7. 


*  Ibid. 
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Figure  6.  Sample  data,  USSPACECOM,  flux  versus  altitude.  May  1990. 
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Figure  7.  Sample  data,  GEODSS,  showing  the  number  of  satellites  seen 
in  an  80.9-hr  period  for  i  >  20  deg. 


•  Impact  data  from  returned  spacecraft  surfaces  were  used  (Refs.  6  and  7).  The  primary 
data  source  was  the  Solar  Maximum  Mission  (Solar  Max)  spacecraft  that  was  in  space 
for  4.15  yr  (1980-84)  at  an  orbital  inclination  of  28.5  deg  and  an  altitude  range  of  500 
to  570  km.  More  than  300  impact  craters  in  the  Solar  Max  thermal  control  louvres, 
representing  ~3  m2  of  returned  space-exposed  surface,  were  examined  (Refs.  6-8). 
Other  impact  data  include  analyses  of  the  windows  from  the  Skylab  (orbital  inclination 
of  50  deg  at  ~400-km  altitude),  Apollo  missions  (moon  missions  and  Skylab 
missions),  and  shuttle  missioiis  (STS-1  to  STS-29,  nominally  at  300-km  altitude). 
Impact  data  must  be  analyzed  to  distinguish  the  namral  micrometeoroid  impacts  from 
the  man-made  (orbital)  debris  impacts  on  the  basis  of  analysis  of  the  residue  in  the 
impact  craters.  It  must  then  be  analyzed  to  determine  the  mass  of  the  impacting 
particles  on  the  basis  of  this  distinction.  Having  made  these  determinations,  the 
nominal  impact  velocities  may  be  estimated.  The  impact  data  were  used  primarily  to 
establish  the  small  particle-size  distribution  for  diameters  of  the  order  of  0.05  cm  and 
less.  Sample  data  are  shown  in  Figure  8. 

•  Recent  radar  measurements  from  Aricebo  and  Goldstone  (Refs.  8  and  9)  are  used  for 
the  size  ranges  of  0.5  to  2  cm  at  altitudes  of  200  to  1000  km  (Aricebo),  and  0.2  to  0.5 
cm  at  a  nominal  altitude  of  575  km  (Goldstone).  These  data  substantiated  the  particle 
size  distribution  for  these  range  of  particles.  Aricebo  data  represent  18  hr  of  viewing 
while  Goldstone  data  represent  14.5  hr  of  viewing.  Sample  data  are  shown  in 
Figure  9. 

•  The  projected  growth  of  mass  in  orbit  for  both  large  and  small  particles  is  derived 
primarily  from  analyses  conducted  using  computer  models  (Refs.  1  and  10),  such  as 
the  NASA  EVOLVE  program  as  well  as  an  analysis  of  the  historical  data.  The  analyses 
examine  both  the  projected  effects  of  various  traffic  models  (primarily  large-particle 
growth)  and  the  effects  the  satellite  breakups  and  random  collisions  (primarily  small- 
particle  growth). 

•  Air  Force  Maui  Optical  Station/Maui  Optical  Tracking  and  Identification  Facility 
(AMOS/MOTIF)  infrared  telescope  data  provided  orbital-debris  particle  albedo  (Refs.  1 
and  5)  used  to  calculate  particle  sizes  from  optical  measurements. 

Comparisons  of  these  data  sets  to  the  NASA  90  model  flux  predictions  were  shown  in  Figures  2 
and  3. 
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Except  for  the  USSPACECOM  and  Solar  Max  data,  the  measurement  periods  represented  by  these 
data  sources  are  extremely  short  The  total  viewing  time  represented  by  the  MTT/LL  ETS, 
GEODSS,  Aricebo,  and  Goldstone  data  is  133.4  hr  or  approximately  5  1/2  days.  This  observing 
time  is  judged  to  be  an  adequate  sampling  time  to  view  a  statistically  significant  number  of  objects 
in  LEO.  Assuming  that  such  a  sample  was  recorded  by  these  measurements,  an  empirically 
derived  flux  model  can  provide  meaningful  results  in  predicting  debris  impacts  on  the  time  scale  of 
years. 

It  has  been  suggested  through  analysis  of  the  GEODSS  data  (Ref.  5)  that  the  USSPACECOM  data 
may  be  incomplete.  Of  the  total  622  satellites  found  by  GEODSS,  only  255  were  identified  or 
found  in  the  USSPACECOM  catalog.  Conversely,  there  are  objects  in  the  catalog  that  were  not 
seen  by  the  GEODSS  telescopes  (Ref.  5).  On  the  basis  of  these  GEODSS  results,  it  was  estimated 
(Ref.  5)  that  the  USSPACECOM-derived  flux  environment  was  low  by  a  factor  of  2  to  4  for  debris 
sizes  around  10  cm  in  diameter. 

The  Solar  Max  data,  which  represent  more  than  4  )rr  exposure  to  the  space  environment,  are  limited 
to  observations  derived  at  only  one  orbital  inclination  and  one  mean  altitude  of  500  km.  Other 
possible  data,  such  as  those  from  the  Infrared  Astronomical  Satellite  (IRAS)  satellite  observations 
(Ref.  11),  were  not  useful  in  the  model  development.  The  IRAS  sensors  and  data  processing 
systems  were  designed  to  filter  out  transient  events,  such  as  orbital  debris  passing  in  view  of  these 
sensors,  especially  at  close  proximities  to  the  IRAS  satellite. 

In  summary,  the  data  used  to  characterize  the  orbital-debris  environment  are  at  best  incomplete  for 
all  LEOs  and  debris  sizes.  Characterization  of  the  debris  environment  can  be  improved  and  data 
uncertainties  reduced  through  a  program  dedicated  to  measure  and  continuously  monitor  the 
complete  orbital-debris  environment  for  all  inclinations,  altitudes,  and  debris  sizes. 

4.2  ASSUMPTIONS 

Certain  key  assumptions  regarding  the  appUcation  of  the  data  base  discussed  in  Subsection  4.1 
were  made  and  led  to  the  development  of  the  NASA  orbital-debris  model  (Refs.  1  and  8).  The  key 
assumptions  to  each  data  set  and  assessment  of  these  assumptions  are  as  follows: 
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•  The  USSPACECOM  data  are  assumed  to  be  complete  to  a  lower  limiting  size  of  10  cm 
at  altitudes  <1000  km,  and  the  distributions  are  assumed  to  be  symmetric  within  each 
latitude  observed.  The  variation  of  flux  with  altitude,  solar  activity,  orbital  inclination, 
and  the  velocity  and  direction  distribution  is  deduced  from  these  data.  The  growth  of 
the  large  particle  distribution  is  recorded  over  an  1 1-yr  solar  cycle  of  the  tracked 
population  of  approximately  5  percent  per  year. 

Assessment:  As  mentioned  in  the  previous  section,  the  USSPACECOM  data  may  not 
be  complete  and  may  be  missing  two  to  four  times  the  actual  number  of  objects  (~10  cm 
diameter)  in  orbit  compared  to  optical  measurements.  This  contention  has  not  been 
settled,  conclusively.  The  diameter  of  the  objects  in  orbit  is  derived  from  the  radar 
cross-section  measurements  (Ref.  12).  These  derived  sizes  depend  on  the  object's 
size,  configuration,  and  orientation  as  well  as  the  radar  system's  characteristics  (Ref. 
13).  The  limiting  diameter  of  10  cm  has  never  been  officially  stated  by 
USSPACECOM  as  the  lower  observing  limit  of  their  systems;  however, 
USSPACECOM  does  refer  to  it  in  presentations  (Ref.  3).  The  assumption  that  the 
USSPACECOM  data  are  complete  down  to  10  cm  must  be  verified.  The  derived 
changes  in  the  flux  environment  with  solar  activity,  etc.,  may  be  valid  even  though  the 
absolute  flux  measurements  may  be  in  error. 

•  The  MIT/  LL  ETS,  during  average-seeing  conditions  (defined  by  observations  of 
known  magnitude  stars),  recorded  twice  the  catalogued  population  ^10-cm  diameter) 
for  >5-cm  particles  and,  during  excellent-seeing  conditions,  recorded  five  times  the 
catalogued  population  for  >2-cm-dia  particles. 

Assessment:  The  size  of  the  observed  object  is  determined  from  its  measured 
brightness  and  assumed  albedo.  Both  depend  greatly  upon  the  object's  shape, 
orientation,  and  surface  properties,  which  are  unknown  for  most  orbital  debris. 
Comparisons  of  GEODSS-derived  diameters  and  radar  cross-section-derived  diameters 
have  resulted  in  an  estimated  average  albedo  for  orbital  debris  of  0.08  with  nearly  half 
of  the  objects  compared  falling  in  the  albedo  range  of  0.05  to  0.20  (Ref.  5).  This  gives 
a  wide  range  of  uncertainty  to  the  actual  size  of  the  object.  To  date,  the  size 
characteristics  for  this  population  have  not  been  verified  with  radar  measurements  for 
fragments  smaller  than  ~8  cm. 
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•  The  measured  impacts  on  the  Solar  Max,  which  included  both  the  meteoroid  and  man¬ 
made  objects,  were  used  as  the  basis  for  calculating  the  orbital-debris  flux.  The  orbital- 
debris  flux  environment  is  1000  times  larger  than  the  meteoroid  flux  for  debris 
diameters  of  >1  ^im,  and  20  percent  of  the  meteoroid  flux  for  debris  diameters  of 
>0.05  cm. 

Assessment:  The  Solar  Max  panels  analyzed  were  exposed  to  the  space  environment 
for  more  than  4  yr  and  should  represent  a  good  sample  of  the  actual  emdronment  at  the 
Solar  Max's  altitude  and  orbital  inclination.  The  distinction  between  meteoroids  and 
man-  made  debris  can  be  established  through  the  analysis  of  the  impacting  particle's 
residue  (Ref.  6).  An  uncertainty  arises  in  determining  the  size  of  the  impacting  object 
without  knowing  the  direction  and  velocity  of  impact  (Normally  the  impacting  object 
is  assumed  to  be  spherical  for  lack  of  any  other  information.) 

•  The  orbital-debris  flux  for  debris  diameters  between  0.05  and  2  cm  is  obtained  by  linear 
interpolation  on  a  log-log  plot  of  flux  versus  diameter. 

AsscssiPgnt*  In  lieu  of  any  other  information  to  the  contrary,  a  straight  line  is  the  most 
conservative  means  of  interpolating,  until  proven  otherwise. 

•  Initially,  the  growth  rate  of  small  particles  was  assumed  to  be  10  percent  per  year  (Ref. 
1),  however,  recent  models  assuming  one  breakup  per  year  suggest  this  growth  rate 
should  be  reduced  to  2  percent  per  year  (Ref.  2).  It  is  assumed  the  major  source  of 
small  debris  is  from  satellite  breakups. 

Assgsgmgnt:  As  will  be  discussed  in  Subsection  4.3.4,  the  historic  record  indicates 
one  breakup  per  year  is  a  reasonable  value.  Modeling  the  resulting  debris  size 
distribution  from  breakups  is  more  difficult;  therefore,  the  growth  rate  of  2  percent 
must  be  validated.  As  will  be  discussed  later,  the  estimated  growth  rates  are  the  most 
controversial  part  of  the  flux  model.  The  assumption  that  the  majority  of  small  debris  is 
from  breakups  is  reasonable  because  the  only  source  for  the  majority  of  the  small  debris 
must  result  from  the  breakup  of  larger  bodies. 

•  The  ratio  of  large-debris  particles  to  small-debris  particles  is  constant  with  altimde,  and 
both  change  at  the  same  rate. 
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Assessment:  Like  the  issue  of  inteipolation  discussed  previously,  there  is  not  enough 
information  on  the  debris  environment,  especially  for  the  small  particles,  to  assume 
otherwise.  Until  more  information  becomes  available,  this  assumption  cannot  be 
verified. 

In  summary,  information  on  the  actual  debris  environment  is  scant,  and  the  ultimate  validity  of  the 
assumptions  used  to  develop  the  model  will  be  verified  only  with  improved  observations.  The 
authors  of  the  NASA  model  are  aware  of  these  limitations  (Refs.  1  and  2)  and  understand  that  the 
model  may  evolve  as  measurements  improve  (as  has  already  happened  with  the  NASA  89  to  the 
NASA  90  model). 

4.3  DEBRIS  FLUX  MODEL  DEVELOPMENT 

With  the  data  base  defined  and  the  assumptions  presented  regarding  the  use  of  this  data  base,  the 
development  of  the  NASA  89/90  model  can  now  be  presented.  Each  functional  relationship  in  this 
model  is  examined.  Specifically,  the  functions  <I>  (h,S),  'F(i),  Fi(d)  and  F2(d),  gi(t)  and  g2(t), 
and  H(d)  are  discussed.  Note  that  the  orbital-debris  model  was  assessed  with  respect  to  the 
parameters  that  were  deemed  the  most  important  (Refs.  1  and  2).  Those  parameters  are  altitude, 
solar  activity,  inclination,  debris  size,  and  time. 

4.3.1  Solar/Atmospheric  Effects.  <I)(h.Sl 

The  function  <I>(h,S),  where 

a>(h,S)  =  <I)i(h,S)  /  [<bi  (h,  S)  +  1]  (NASA  89/90)  (15) 


and 


<I>i(h,S)  =  10(‘>/200-s/14o-i.5)  (NASA  89/90)  (16) 

accounts  for  the  influence  of  the  upper  atmosphere  in  "cleansing"  or  removing  the  orbital  debris 
because  of  atmospheric  drag  effects.  It  also  demonstrates  the  influence  of  the  solar  activity  in 
increasing  or  decreasing  this  "cleansing"  rate.  The  function  ^>(h,S)  remains  unchanged  between 
the  NASA  89  (Eq.  1)  and  NASA  90  (Eq.  6)  models.  The  behavior  of  O  with  altitude,  h,  and  the 
10.7-cm  wavelength  measure  of  the  solar  spectrum,  S,  (defined  for  the  previous  year,  i.e.,  t  - 1, 
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and  shown  to  be  an  accurate  measure  of  solar  activity  along  with  sunspot  counts)  is  shown  in 
Figure  10.  Note  that  this  function  is  normalized  such  that  its  maximum  value  is  one  and  modifies 
the  flux  magmtude  as  altitude  and  solar  activity  varies.  At  any  given  altitude,  when  the  solar 
activity  increases  (increasing  S),  <J>  decreases,  reducing  the  flux  environment .  As  altitude 
increases  for  a  given  value  of  S,  the  effects  of  the  atmosphere  decrease.  The  maximum  changes 
with  respect  to  S  appear  at  =500  km.  Beyond  some  altitude  limit,  nominally  1000  km,  the  effects 
of  the  atmosphere  are  negligible.  The  4>  appears  to  behave  properly. 


Figure  10.  Solar/atmospheric  effects,  ^>(h,S). 

The  extent  of  upper  atmosphere  where  orbiting  objects  are  affected  by  solar  activity  is  nominally  up 
to  800  km  (Refs.  14  and  15).  The  effect  of  the  atmospheric  drag  is  to  reduce  the  energy  of  the 
orbiting  body  by  first  making  the  orbit  circular.  Once  the  orbit  is  circular,  it  then  begins  to  decay 
leading  evenmally  to  the  object  re-entering  the  atmosphere.  It  can  be  shown  that  the  lifetime  of 
objects  in  circular  orbits  is  inversely  proportional  to  the  atmospheric  density  (Ref.  15).  One 
measure  of  solar  activity  is  characterized  by  the  10.7-cm  wavelength  of  the  solar  spectrum 
(Refs.  15  and  16)  averaged  over  a  typical  13-mo  period  of  the  11-yr  solar  cycle  (Ref.  16). 
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Nominal  values  of  S  (10^  Jy)  during  a  solar  cycle  range  from  a  value  of  70  x  10^  Jy  to  200  x  10^ 
Jy;  however,  values  exceeding  200  x  10^  Jy  and  approaching  250  x  10^  Jy  have  been  observed 
(Refs.  15  and  16).  The  implication  and  result  of  wide  variations  in  the  maximum  value  of  S  during 
a  solar  cycle  make  long-term  debris  predictions  difficult  (Refs.  16  and  17). 

The  formulation  of  «J>  is  shown  in  Figure  11.  The  <!>  is  a  composite  function  made  up  of  <I>i(h,S) 
and  ^>2  (h,S).  The  function  Oi  models  the  effects  through  the  altitudes  where  atmospheric  effects 
are  important.  The  ^  models  the  debris  flux  at  altitudes  beyond  the  upper  atmosphere.  From 
Figure  11,^  was  assumed  to  predict  a  constant  flux  at  these  altitudes.  The  function  was 
determined  by  a  curve  fit  to  the  USSPACECOM  data  through  the  lower  altitudes  where  solar 
influences  on  the  debris  environment  are  important  The  curve  fit  to  the  data  shows  an  increase  in 
the  environment  by  an  order  of  magnitude  for  every  200-km  altitude  increment;  thus  the  factor, 
h/200.  The  nominal  level  of  solar  activity  was  picked  to  be  140  x  Ifl^  Jy  and  the  function  given 
by  Equation  16  resulted  from  these  curve  fits.  Because  the  actual  magnitude  of  the  debris  flux 
environment  is  contained  in  Fi  and  Fj  (times  the  growth),  the  function  4>  2  was  set  equal  to  one  to 
normalize  ^>.  The  function  4>  modifies  the  predicted  flux  environment  (magnitude),  F,  produced 
from  the  sum  of  the  small-  and  large-particle  fluxes,  Figi  +  F2g2- 


Figure  11.  Construction  of  <I>(h,S). 

4.3.2  Inclination  Distribution.  ^0) 

Although  the  orbital-debris  flux  environment  is  assumed  to  be  symmetric  in  latitude  (evenly 
distributed  around  the  earth),  observations  show  that  it  varies  as  a  function  of  orbital  inclination 
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(Ref.  18).  Typical  data  showing  these  variations  are  depicted  in  Figure  12.  The  function  'F(i) 
adjusts  the  predicted  flux  environment  accordingly  as  a  function  of  the  orbital  inclination.  The  'P  is 
shown  in  Figure  13  (the  corresponding  values  were  shown  in  Table  1).  The  function  'P,  in  curve 
fitting  the  data,  is  normalized  for  a  uniform  random  distribution.  The  peaks  (polar  orbits)  represent 
inclinations  where  there  is  more  traffic  and  where  several  on-orbit  breakups  have  occurred.  The 
overall  validity  of  'P  rests  with  the  validity  of  the  measured  concentrations  of  debris.  The  most 
serious  concern  with  'P  is  that  it  does  not  account  for  different  altitude  variations  with  i. 
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Figure  12.  Inclination  distribution  of  USSPACECOM-tracked  obiects 
August  1988  -  all  altitudes. 


versus  orbital  inclination,  i  (deg). 
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4.3.3  Particle  Rux.  Fi  fd)  and  F7  (d).  H  fd) 


The  prediction  of  the  orbital-debris  flux  environment  is  given  by  Fi  for  small  particles  and  F2  for 
large  particles.  The  nominal  distinction  between  small  and  large  particles  is  in  the  range  of  1-  to 
10-cm  diameter.  Figure  14  shows  how  Fi  and  F2  vary  with  diameter.  As  seen  in  this  figure,  Fi 
dominates  for  the  smaller  particles  and  F2  for  the  larger  particles.  They  intersect  at  ~3.2-cm 
diameter,  which  is  in  the  range  where  the  distinction  between  large  and  small  particles  is  made. 


Figure  14.  Small-particle  flux,  Fi(d),  and  large-particle  flux, 
F2(d),  versus  diameter,  d(cm). 


The  small-particle  flux  distribution,  Fi,  was  derived  using  data  from  two  sources:  impact  data 
from  returned  spacecraft  surfaces  and  optical  data  (ground  observations).  The  impact  data  (Refs. 
6-8)  consist  primarily  of  analyses  of  the  returned  surfaces  (thermal  control  louvres)  from  the  Solar 
Max  satellite,  which  operated  at  a  nominal  altitude  of  5(X)  km  and  orbital  inclination  of  28.5  deg. 

In  the  analysis  to  determine  the  orbital-debris  environment,  it  was  assumed  the  meteoroid  flux  that 
impacted  the  Solar  Max  was  a  known  percentage  of  the  total  flux.  The  orbital-debris  environment 
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is  determined  from  the  ratio  of  meteoroid  to  orbital-debris  craters/impacts  as  presented  in 
Subsection  4.2.  The  measured  impact  data  resulted  in  defining  the  orbital-debris  flux  for  diameters 
nominally  up  to  0.05  cm.  The  optical  data  from  the  MIT/LL  ETS  (Ref.  4)  was  analyzed  by 
NASA* .  The  analysis  derived  from  these  data  form  the  "anchor"  for  the  curve  fit  at  the  larger 
diameters  of  2  and  5  cm.  Without  any  other  data,  the  curve  fit  between  0.05  and  2  cm  is  assumed 
to  be  linear  Oog-log  relationship  of  flux  versus  diameter).  It  must  be  noted  that  ~8  cm  is  the 
smallest  debris  size  whose  population  can  be  reasonably  confirmed  with  radar  measurements. 

The  USSPACECOM  data  were  used  to  determine  the  large-particle  flux  distribution,  F2.  After  the 
formulation  of  the  NASA  89  model,  analysis  of  GEODSS  data  led  NASA  to  an  increase  of  the 
predicted  environment  in  the  10-cm-dia  range  (Ref.  6).  As  discussed  in  Subsection  4.2,  analysis 
of  the  GEODSS  optical  data  in  the  8-  to  30-cm-dia  range  appears  to  have  revealed  two  to  four  times 
the  number  of  objects  seen  by  the  USSPACECOM-catalogued  radar  data  (Ref.  5).  The  "Henize" 
function,  H,  represents  a  model  adjustment  that  increases  the  environment  in  this  range.  Figure  15 
shows  the  behavior  of  H  versus  diameter.  It  is  a  correction,  peaking  at  a  value  of  3.16  for  a 
diameter  of  6.025  cm,  and  increasing  the  predicted  environment  for  the  1-  to  100-cm  dia  debris 
particles.  Although  the  data  only  support  corrections  above  8-  to  10-cm  dia  particles,  the  H 
correction  is  applied  to  smaller  diameters  and  represents  an  extrapolation  to  the  smaller  sizes.  The 
Aricebo  and  Goldstone  radar  data  (Ref.  9)  confirm  the  NASA  89  model  predictions  in  the  range  of 
0.2-  to  0.5-cm  diameter  (Goldstone)  and  0.5-  to  2-cm  diameter  (Aricebo),  within  the  measurement 
uncertainties. 


Figure  15.  Flux  correction,  H(d),  versus  diameter,  d(cm)  for  NASA  90  model. 


*  Internal  NASA/JSC  Memo  by  Loretta  Weiss,  "Reanalysis  of  the  MIT  ETS  Telescope  Data,"  1985. 
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Comparisons  of  the  NASA  90  model  predictions  to  data  from  Solar  Max,  Goldstone,  Aricebo,  and 
GEODSS  were  shown  in  Figure  2  (model  flux  versus  diameter,  h  =  500  km)  and  to 
USSPACECOM  data  in  Figure  3  (model  flux  versus  altitude,  d  ^  10  cm).  Figure  2  shows  the  fit 
to  the  small-particle  data  while  Figure  3  shows  the  fit  to  the  large-particle  or  USSPACECOM  data. 
The  curve  fit  to  the  small-particle  data  in  Figure  2  shows  good  agreement  since  these  data  are  the 
basis  for  the  curve  fit  Figure  3,  shows  a  good  fit  to  the  USSPACECOM  data  at  the  lower 
altitudes,  but  illustrates  the  effect  of  the  low-order  curve  fit  and/or  lack  of  confidence  in  the  data  at 
higher  aitimdes.  The  obvious  peaks  and  valleys  in  the  data  at  the  higher  altitudes,  h  >  800  km,  are 
not  represented  by  the  curve  fit 

The  fluctuation  in  the  data  at  -800, 1000,  and  1500  km  are  thought  to  be  real.  They  are  due  to 
known  breakups  at  these  altitudes.  Approximately  150  such  events  (either  intentional  or 
unintentional)  have  been  recorded  to  date  at  various  altitudes  (Refs.  19  and  20).  Some  of  the  more 
severe  breakups  are  shown  in  Table  2  (Refs.  19  and  20). 

4.3.4  Debris  Growth.  gi(t)  and  goft) 

The  most  controversial  part  of  the  debris  flux  model  concerns  the  predicted  growth  of  the  space- 
debris  envirorunents  with  time.  The  controversy  concerns  the  magnitude  and  whether  the  growth 
rate  is  linear  or  compounded.  Although  the  growth  rate  may  be  bracketed  on  the  basis  of  historical 
data  and  model  simulations,  there  simply  is  not  enough  data  to  confirm  or  quantify  the  rate.  In  the 
revision  from  the  NASA  89  to  the  NASA  90  model,  a  major  change  is  seen  in  the  form  of  the 
growth  functions,  gi(t)  and  g2(t),  and  the  growth  rates,  p  and  q.  The  differences  are: 

Function 
gi(t) 
g2(t) 

P 

q 


NASA  89 

(1  +  2p)  (H985) 
(1  +  P)  (t-1985) 
-5% 

N/A 


NASA  90 

(1  +  q)  (t-1988) 

l+p(t-1988) 

-5% 

-  2%  until  2010, 
then  4%  thereafter 
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Table  2.  Fragmentations/breakups. 


Object 

(Common  Name) 

Breakup 

Date 

Approximate 
Altitude  (km) 

Catalogued  Objects  (d  >  10  cm) 
(Approximate) 

Initially  Still  in  Orbit  (1988) 

Titan  3C-4 

1965 

740 

470 

90 

Ariane  3rd  Stage 

1986 

820 

465 

460 

Thor-Agena-D 

1970 

1080 

345 

295 

Cosmos  1275 

1981 

980 

300 

290 

Ablestar  Rocket 

1961 

950 

270 

210 

Thor-Agena-D 

1969 

920 

265 

140 

Solwind 

1985 

530 

250 

120 

Cosmos  844 

1976 

210 

250 

0 

Delta  2nd  Stage 

1975 

725 

230 

90 

Delta  2nd  Stage 

1976 

750 

200 

50 

Cosmos  544 

1973 

310 

200 

0 

Delta  2nd  Stage 

1981 

900 

195 

175 

Delta  2nd  Stage 

1973 

1500 

185 

170 

Cosmos  1813 

1987 

350-410 

185 

50 

Cosmos  57 

1965 

160-700 

165 

0 

Delta  2nd  Stage 

1977 

1450 

160 

100 

Cosmos-Rocket 

1965 

1640 

150 

25 

Delta  2nd  Stage 

1975 

1460 

140 

135 

Delta  2nd  Stage 

1977 

1510 

140 

135 

Cosmos  252 

.  1968 

530 

130 

60 

The  small-particle  growth  rate  has  been  decoupled  from  the  large-particle  growth  rate  through  the 
introduction  of  the  coefficient,  q,  in  the  NASA  90  model.  Previously  it  was  assumed  the  small- 
particle  growth  rate  was  twice  that  for  large  particles  (Ref.  1).  Although  the  growth  of  small 
particles  is  still  compounded  with  time,  the  rate  (2  percent  for  the  NASA  90  versus  10  percent  for 
the  NASA  89)  is  significantly  less  in  the  NASA  90  model.  Next,  the  large-particle  growth  has 
been  changed  from  a  compounded  growth  to  a  linear  growth  with  time.  The  assumed  growth  rate, 
p  =  5  percent,  remains  the  same.  The  overall  result  of  these  changes  is  to  reduce  the  projected 
growth  rate  with  time  as  shown  in  Figure  16,  especially  for  the  small  particles. 
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Figure  16.  Small-particle  growth,  g  i  (t),  and  large-particle  growth,  g  2  (t),  NASA  89 
and  NASA  90  models  versus  time,  t  (yr). 

The  basis  for  these  projections  is  historical  data  such  as  depicted  in  Table  2  and  model  simulations 
to  estimate  the  amount  of  small  (untrackable)  debris  produced.  The  major  source  of  smaU-particle 
debris  growth  is  assumed  to  result  from  satellite  breakups.  Over  the  past  decade,  there  have  been 
37  breakups  in  LEO  producing  10  or  more  trackable  fragments.*  Of  these  37  breakups,  11  were 
thought  to  be  unintentional,  and  of  these  11, 8  were  felt  to  be  accidental  explosions.  The 
remaining  three  were  possible  collisions.  (The  Cosmos  1275  is  considered  to  be  a  strong  collision 
candidate.)  Assuming  in  the  future  that  intentional  breakups  are  reduced  by  improved  operating 
procedures  and  heightened  awareness  by  users  not  to  create  more  long-term  debris,  one  can 
assume  (based  on  the  historical  data)  the  rate  of  one  unintentional  breakup  per  year  will  continue. 
On  the  basis  of  this  breakup  rate  and  the  NASA  EVOLVE  model  predictions,  a  2-percent  growth 
rate  was  estimated.  This  is  felt  to  be  a  valid  approximation  through  the  year  2010.*  Small-debris 
growth  is  still  assumed  to  be  compounded  based  on  simulations.  If  current  practices  continue,  q  is 
projected  to  increase  to  4  percent  on  or  about  2010  (Ref.  2). 


* 


Private  meeting s/communications  with  NASA/JSC,  DJ.  Kessler,  1990-1991. 
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The  growth  of  large  particles  is  found  in  historical  records  and  is  related  primarily  to  launch  rates 
and  objects  in  orbit  These  data,  recorded  by  USSPACECOM,  were  used  to  determine  and  to 
project  the  future  growth.  Data  are  averaged  over  solar  cycles.  During  the  1966-1977  solar  cycle, 
the  average  increase  in  the  catalogued  population  was  300  objects  per  year,  representing 
approximately  a  5-percent  increase  per  year.  This  is  the  basis  for  the  value  of  p.  The  fumre 
activity  and  manner  in  which  the  large-particle  debris  population  will  increase  cannot  be  deduced 
from  this  historical  record.  It  is  clear  from  Rgure  17  that  the  data  do  not  justify  either  linear  and 
compounded  growth.  Traffic  models  for  a  variety  of  conditions  were  examined  to  aid  in  predicting 
and  bracketing  the  possible  future  growth  trends  (Refs.  1  and  2),  from  constrained  to  high  world- 
traffic  models.  The  linear  growth  was  chosen  based  on  a  continued  350  launches  per  year  and 
assuming  a  1988  catalogued  population  of  1000. 


Figure  17.  Historical  large-particle  growth. 
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4.4  DEBRIS  IMPACT/DTRECTTON  DTSTRTBTJTION  AND  PARTICLE  MASS 


4.4.1  Colli.sion-Velocitv  Impact  Distribution 

The  expression  for  the  collision-velocity  impact  distribution,  f(V),  as  shown  in  Equation  8,  gives 
the  number  of  impacts  for  orbital  debris  having  velocities  between  V  and  V  -t-  dV  relative  to  a 
spacecraft  with  orbital  inclination,  i.  The  normalized  collision- velocity  impact  distribution  is  given 
by  Equation  10.  The  collision-velocity  impact  distribution  was  determined  from  USSPACECOM 
data  for  altitudes  between  500  and  1000  km  (Refs.  1  and  2).  Examples  of  those  data  are  shown  in 
Figure  18  (Ref.  10).  The  data  were  examined  for  variations  with  parameters  such  as  time,  altitude, 
or  inclination.  The  only  significant  variation  that  could  be  distinguished  was  the  variation  with 
respect  to  the  orbital  inclination  (Refs.  1  and  2).  Variations  both  in  time  and  altitude  appeared  to 
look  like  random  variations  and  were  regarded  as  second-order  effects*.  This  led  NASA  to 
formulate  f  (V),  given  in  Equation  8. 

The  function  f  (V)  is  a  double  Guassian  fit  to  the  data  in  Figure  18.  Although  there  may  be  some 
size  dependency,  data  are  insufficient  to  distinguish  this  effect*.  The  function  f  (V)  is  a 
complicated  expression  without  any  physical  significance  and  whose  empirically  determined 
constants  vary  with  inclination  in  a  seemingly  unrelated  fashion.  Finally,  it  is  not  an  expression 
that  is  easily  used. 


Figure  18.  Velocity  distribution  data  (h  =  500  km,  i  =  28.5  deg). 


*  Ibid  p.  30. 
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4.4.2  Impact  Direction  Distribution 


The  distribution  of  debris  impact  direction  (Eq.  1 1)  in  a  specified  orbital  inclination  is  symmetric 
and  does  not  account  for  out-of-plane  impacts.  On  the  basis  of  measured  data,  the  average  velocity 
for  the  debris  was  found  to  be  7.7  km/s;  this  represents  an  ideal  circular  orbital  altitude  of  345  km 
while  the  data  were  averaged  over  500  to  1000  km. 

NASA  has  deduced  that  there  are  indeed  debris  with  elliptical  orbits.  This  results  in  an  average 
velocity  of  7.7  km/s,  which  is  higher  than  expected  from  debris  only  in  circular  orbits.  A 
rearward-facing  surface  (surface  normal  opposite  the  direction  of  motion)  of  an  object  in  a  circular 
orbit  would  not  be  impacted  by  debris  unless  there  were  debris  in  elliptical  orbits.  Returned 
surfaces,  including  the  Long  Duration  Exposure  Facility,  show  impacts  on  the  rear  surfaces,  thus 
substantiating  that  there  exists  debris  in  elliptical  orbits.  There  are  definitely  out-of-plane  impacts 
as  well  as  impacts  on  rear  surfaces,  thus  the  expression  (Eq.  1 1)  represents  only  an  approximation 
of  the  actual  in-plane  debris  environment 

4.4.3  Orbital-Debris  Particle  Density 

The  expression  for  the  orbital-debris  particle  density  assumes  that  the  particles  are  spherical  in 
shape  (Eq.  12).  There  is  little  information  on  actual  debris  shapes.  There  does  exist  a  data  base 
for  particle  sizes  >0.5  cm  which  consists  of  studies  of  orbital  decay  and  fragmentation  experiments 
involving  typical  spacecraft  structures  (Refs.  18, 21,  and  22).  These  data,  as  described  in 
Reference  21,  have  led  to  the  expression  (Eq.  12)  for  particle  density  (d  ^  0.5  cm),  which  is 
assumed  to  be  accurate  within  a  factor  of  0.5  to  2  (Ref.  2). 

For  sizes  <0.5-cm  diameter,  the  data  are  virtually  nonexistent  It  is  known  that  orbital  debris 
consists  of  a  variety  of  material  types,  including  epoxy-glass  (p  ~  1.8  g/cm  3),  aluminum 
(p  ~  2.8  g/cm  3),  copper  (p  ~  8.9  g/cm  3),  and  steel  (p  -  8  g/cm  3).  Estimates  could  be  made  of 
the  approximate  composition  of  the  debris  (by  volume  fraction)  because  of  fragmentation, 
collision,  and  operational  procedures  (including  engine  operations);  however,  definitive  studies 
have  not  been  performed  (Ref.  2).  The  value  of  4.7  g/cm  3  for  d  <  0.5  cm  results  from 
Equation  12  at  d  =  0.5  cm. 

The  large  uncertainty  in  particle  mass  has  implications  for  impact  and  shielding  studies.  For  impact 
and  shielding  studies,  it  is  important  to  know  not  only  the  energy  or  momentum  of  the  impacting 
particle,  but  also  the  mass  per  unit  area  of  the  impactor*.  Mass  is  extremely  important  to  satellite 

Private  meetings/communications  with  PLyWSSD,  Dr.  Firooz  AUahdadi,  1990-1991. 
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designers  (who  want  to  minimize  mass),  and  unnecessary  spacecraft  mass  due  to  unneeded 
shielding  results  in  increase  system  costs.  In  using  an  expression  for  debris  mass  that  may  be  in 
error  by  a  factor  of  2,  a  shield  design  may  have  either  too  little  shielding  or  too  much.  Either 
extreme  represents  a  dilemma  to  the  satellite  designer.  The  present  expression  for  density  along 
with  the  collision-velocity  impact  distribution  and  direction  predictions  is  felt  to  represent  extremely 
gross  approximations  of  the  environment. 

4.5  CIRCULAR  VERSUS  RTJ.IPTICAL  ORBITS 

In  the  NASA  89/90  models,  both  the  orbital  debris  and  the  spacecraft  of  interest  are  assumed  to  be 
in  circular  orbits  for  the  purpose  of  computing  the  direction  of  impact,  a  (Eq.  1 1).  Because  V  is 
always  positive.  Equation  1 1  allows  for  impact  directions  between  only  ±  90  deg  (0  deg  being  the 
spacecraft  direction  of  motion).  The  formulation  for  impact  direction  does  not  consider  out-of¬ 
plane  (parallel  to  Earth's  surface)  impacts.  Equation  1 1  is  plotted  in  Figure  19  for  the  range  of 
debris-impact  velocities,  0  to  15.4  km/s.  The  maximum  impact  velocity  occurs  at  a  =  0  deg.  The 
impact  velocity  is  zero  at  a  >  ±  90  deg.  Solving  Equation  1 1  for  V  and  combining  with 
Equation  8,  the  expected  number  of  impacts  from  any  direction  can  be  calculated.  This  results  in 
the  "butterfly"  diagram  shown  in  Figure  20.  ITus  figure  emphasizes  the  directionality  of  impacts 
because  of  the  orbital-debris  flux  environment  (Refs.  22  and  23). 


0° 


Figure  19.  NASA  90  model  debris-impact  velocity  envelope,  versus  direction  of 
impact,  a,  where  the  spacecraft  is  in  the  center  of  polar  coordinates. 
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Figure  20.  Debris-impact  envelope,  V  •  f  '(V),  NASA  90  model  (i  =  28.5  deg). 

4.6  DATA/MQDEL  UNCERTATNTTFS 

The  orbital-debris  engineering  model  developed  by  NASA  is  based  on  a  curve  fit  to  debris 
measurement  data  (Refs.  1  and  2).  The  accuracy  of  this  type  of  empirical  model  is  related  to  the 
fidelity  of  the  curve  fit  and  the  representativeness  of  the  data  sample  used.  The  uncertainties  due  to 
curve  fits,  measurement  resolution  and  sampling  techniques  have  been  estimated  by  NASA 
(Ref.  2)  and  are  shown  in  Table  3.  Values  in.this  table  will  form  the  basis  for  the  uncertainty 
analysis  to  be  presented  in  Section  7.0. 


The  uncertainties  given  in  this  table  represent  the  90-percent  confidence  estimated  values*  based  on 
the  consensus  of  NASA  investigators.  The  first  parameter,  the  uncertainty  in  the  flux 
measurements,  reflects  statistical  and  measurement  uncertainties  in  the  data  set  for  the  various  sizes 
of  debris.  The  uncertainty  in  the  altitude  distribution  contributes  to  the  total  uncertainty.  The 
uncertainty  is  based  on  the  difficulty  in  including  flux  for  debris  in  highly  elliptical  orbits.  No 
altitude  distribution  uncertainties  are  presented  in  the  1-  to  10-cm-dia  size  range.  The  uncertainty  in 
the  collision-velocity  impact  distribution,  f(V),  was  estimated  for  the  fraction  of  debris  with  impact 
velocities  <5  km/s.  The  uncertainties  for  debris  density  and  shape  reflect  the  extreme  lack  of 
knowledge  concerning  small  particles,  d  <  1  cm.  The  variations  presented  for  growth  rates  and 
level  of  solar  activity  represent  extreme  values. 


* 


Private  meetings/communications  with  NASA/JSA,  D.J.  Kessler,  1990-1991. 
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Table  3.  Estimated  orbital-debris  uncertainties,  NASA  90  model,  90-percent  confidence. 


Parameter 

Estimated  Uncertainty 

Flux  Measurements 

d  >  10  cm 

1.5  to  0.5  X  Flux 

0.05  <  d  <  10  cm 

3.0  to  0.33  X  Flux 

d  <  0.05  cm 

2.0  to  0.5  X  Flux 

Altitude  Distribution 

d  >  10  cm 

2.0  to  0.5  X  Flux 

d<  1  cm 

5.0  to  0.2  X  Flux  per  every  200  km  away  from 
h  =  5(X)km 

Velocity  Distribution 

V  <  5  km/s 

0.5  to  3.0  X  f(V)  only  for  fraction  where 
V<5km/s 

Debris  Density 

d  >  1  cm 

2.0  to  0.5  X  mean  density 

d  <  1  cm 

Unknown,  can  estimate  by  estimating 
volume/weight  fractions  or  use  +100%  for 

p  =  2.8  g/cm3 

Debris  Shape 

Unknown 

Orbital-debris  Growth 

P 

4  to  10% 

q 

0to20% 

Solar  Activity 

S 

Use  max  and  min  values  for  nominal  solar 
cycle  where  typically  70  <  S  <  210 

4.7  MODEL  RESTRICTTQNS/LIMITATTON.S 

Initial  examination  of  the  model  development  and  of  the  corresponding  assumptions  resulted  in  the 
documentation  of  the  qualitative  model  limitations.  The  limitation  for  use  of  the  model  to  LEO  was 
noted  by  NASA  (Refs.  1  and  2).  The  use  of  the  NASA  90  model  should  be  limited  in  the  temporal 
extent  of  its  use.  The  time  validity  of  the  model,  based  on  current  projections,  is  assumed  to 
extend  to  the  year  2010. 

General  limitation  is  that  this  model,  developed  primarily  using  historical  data,  reflects  the  way 
space  operations  have  been  conducted.  If  one  drastically  changes  current  operational  procedures  or 
traffic  model  assumptions,  then  the  model  may  no  longer  be  applicable.  Use  of  this  model  must  be 
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consistent  with  the  fundamental  assumptions  concerning  its  development  It  is  a  tool  for  use  by 
design  engineers  in  estabhshing  an  order  of  magnitude  orbital-debris  hazard.  It  is  not  meant  for 
sophisticated  trade  studies  or  examinations  of  "what  ifs"  in  space  architecture  studies. 

Finally,  it  is  clear  that  the  model  does  not  fit  precisely  all  fluctuations  in  the  data,  most  notably 
those  altitudes  and  inclinations  where  known  breakups  have  occurred.  If  one  specifically  wishes 
to  predict  debris  environments  with  higher  accuracy,  then  a  debris  model  with  a  higher  order  fit  to 
the  data  is  needed.  The  model  will  be  required  to  simulate  and  integrate  the  effects  of  breakups  to 

improve  its  accuracy.  Model  capabilities  and  limitations  will  be  further  quantified  in  Sections  6.0 
and  7.0. 


37 


5.0  THE  AF/DoD  ORBITS  OF  INTEREST 


The  AF/DoD  orbits  of  interest  have  been  identified  on  the  basis  of  an  analysis  of  the  historical 
record  and  projected  mission  requirements  (Ref.  24).  From  this  analysis,  it  appears  that  future 
AF/DoD  traffic  will  operate  in  five  orbital  regimes  (Ref.  24)  as  shown  in  Table  4  and  depicted  in 
Figure  21.  These  projections  are  based  on  the  analysis  of  all  AF/DoD  requirements  through  the 
year  2010.  If  there  is  any  preponderance  of  new  traffic,  it  will  be  destined  for  Regimes  D  and  E  at 
relatively  low  altitudes  (h  ^  1200  run  or  h  <  2220  km).  Note  that  the  Space  Station  Freedom  is 
destined  for  operation  in  Regime  C  (460-km  altitude,  28.5-deg  orbital  inclination). 


Table  4.  The  AF/DoD  orbits  of  interest. 


Regime 

Description 

Orbital 

Inclination 

(deg) 

Altitude  (h) 

A 

High  Altitude/ 
Geosynchronous 

0<i<67 

-19330  run  (-35800  km) 
[Synchronous] 

B 

Mid  Altitude/ 
Midinclination 

55  <  i  <  67 

5600nm<h^  11000  run 
(10400  km  <h^  20400  km) 

C 

Low  Altitude/East 

28  <  i  <  32 

h  <  1000  run  (h  <  18540  km) 

D 

Low  Altitude/ 
Midinclination 

60<i<80 

h  <  1000  run  (h<  1850  km) 

E 

Low  Altitude/Polar 

90<i<  100 

h  <  4000  run  (h<  7400  km) 

On  the  basis  of  the  analysis  of  the  preceding  sections,  NASA  90  orbital-debris  prediction  for  the 
AF/DoD  orbits  of  interest  will  be  made  only  for  those  regimes  where  the  results  are  applicable  as 
shown  in  Figure  21:  RegimesC,  D,  E  for  h<  2000  km. 
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Regime  A 


GEO 


Figure  21.  The  AF/DoD  orbits  of  interest. 
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6.0  SENSITIVITY  STUDY 


A  sensitivity  study  of  the  NASA  90  model  is  discussed  in  this  section.  The  sensitivity  or  variation 
of  the  orbital-debris  flux  prediction  is  calculated  by  examining  the  first  partial  derivative  of 
F  (d  ,h  ,i  ,t  ,S)  with  respect  to  the  variables  d,  h,  i,  t,  S,  p,  and  q.  The  derivative  gives  the  relative 
change  of  F  with  respect  to  the  change  in  each  pai^eter.  The  purpose  of  the  sensitivity  study  is  to 
identify  those  variables  that  cause  the  greatest  change  in  F  for  given  changes  in  the  variables.  The 
more  sensitive  parameters  will  cause  the  greatest  change  in  F.  Any  anomalous  behavior  is  also 
examined,  especially  any  singularity  in  the  function,  F.  The  derivative  equations  are  presented  in 
Subsection  6.1  while  the  analyses  are  presented  in  Subsection  6.2.  Subsection  6.3  summarizes  the 
results  of  the  sensitivity  study. 

6.1  SENSmVrrY  parameters/analyttc  fxprfsstons 

The  function  F(d,  h,  i,  t,  S)  is  composed  of  several  distinct,  empirically  derived  functions. 

Because  the  majority  of  functions  that  make  up  F  are  smooth  and  continuous,  it  is  expected  that  the 
first-order  derivatives  also  will  be  smooth  and  continuous.  The  function  'F(i)  is  based  on  tabular 
values,  which  when  plotted  results  in  a  smooth  continuous  curve  over  the  range 
25  deg  ^  i  <;  125  deg.  The  function  gi(t)  is  piecewise  continuous  with  the  discontinuity  occurring 
where  q  changes  value  at  the  year  2010.  The  functional  forms  making  up  F  (d,  h,  i,  t,  S)  are 
shown  in  Figures  10  ^d  13  through  16  for  ^(h,  S),  'P(i),  Fi(d)  and  F2(d),  H(d)  and  gi(t)  and 
g2(t),  respectively. 

The  first-order  partial  derivatives  of  F  with  respect  to  each  variable  are  as  follows: 


—  =(-i? 
3d 


l{log,„d-0.7i)  ^-1 

log  ,od- 0.78]  ^  J 
0.637  /  _j_| 

F.S.  i 

[ _ 6_] 

1  d  +  700J 

1  F2g2l 

d(  0.637)  2 

t  Figi  +  F2g2 

/ 

(17) 


3Fj3y/3i)p^(Ay/Ai)j, 
ai  'F  ~  T' 

where  A'F/Ai  is  determined  by  numerical  finite  differencing. 


(18) 


(19) 


40 


_  [Figi  In  (1  +  q)  +  PF2I  p 

at  Figi  +  F2g2  J 


(20) 


dp 


(t-  1988)F, 

Figi+F2g2. 


F 


dq 


'(t-  1988) 
(1  +  q) 


PiSi 


L  Figi  +  F2g2 


F 


(21) 


(22) 


(23) 


6.2  SENSITIVITY  ANALYSTS  AND  RF..STTT  .T<; 

In  performing  the  sensitivity  analyses,  variations  of  the  parameters  for  two  base  line  cases  were 
made  as  follows: 


CaS£i:  d  =  1  cm,  h  =  500  km,  i  =  47  deg  ('F  (i)  =  1.0),  t  =  1988,  S  =  140,  p  =  0.05, 

and  q  =  0.02 

d  =  1  cm,  h  =  500  km,  i  =  47  deg,  t  =  1995,  S  =  90,  p  =  0.05,  and  q  =  0.02 

There  is  no  particular  significance  in  the  choice  of  these  two  cases  other  than  selecting  a  value  of 
unity  for  debris  diameter  and  the  orbital  inclination.  One  selects  the  model  starting  year,  1988,  and 
some  future  year,  1995,  to  ensure  that  significant  changes  with  respect  to  time  are  identified  (the 
values  of  S  are  representative  of  those  years).  Finally,  the  values  of  p  and  q  represent  the  model 
nominal  values. 


From  these  base  line  values,  variations  with  respect  to  each  derivative  functional  relationship  were 
calculated  for  the  following  ranges: 


•  —  versus  d  for  10'^  <  d  <  10^  cm 

9d 


•  — versus  h  for  1(X)  <  h  <  2(X)0  km 
3h 
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ap 

•  — versus  i  for  25  <  i  <  125  deg 

ai 


ap 

•  —  versus  t  for  1988  or  1995  <  t  <  2010 

at 

ap 

•  —  versus  S  for  70  <  S  <  250 

as 

ap 

•  —  versus  p  for  0  <  p  <  0.20 
ap 

ap 

•  —  versus  q  for  0  <  q  <  0.20 
aq 


6.2.1  ap/ad  versus  d 


The  variation  of  ap/ad  versus  diameter  is  shown  in  Pigure  22  for  Cases  I  and  n  and  Pigures  23a 
and  b  for  Case  I.  Prom  Pigure  23, 


i.e.,  the  flux  becomes  more  sensitive  to  particles  with  smaller  diameters.  Note  that  there  are  16 
orders  of  magnitude  change  in  ap/ad  for  the  range  of  diameters,  0.001  to  100  cm.  The  rapid 
change  of  ap/ad  is  especially  apparent  when  plotted  on  linear  axes  in  Pigure  22b.  The  expression 
lap/adl  =  10^  at  d  =  0.001  cm,  10-  at  d  =  1  cm,  and  10-  ^  at  d  =  100  cm.  Note  that  lap/adl  <1  for 
diameters  >0.05  cm.  There  is  a  slight  change  in  slope  of  ap/ad  versus  diameter  at  d  =  5.0  cm. 
Por  diameters  <0.1  cm,  Pj  (and  thus  its  derivative)  dominates  while  for  diameters  >10  cm,  P2 
dominates;  thus,  the  reason  for  the  change  in  slope.  Note,  in  the  limits  (for  the  absolute  values), 

- >  00  as  d  0 

ad 


This  means  that  the  flux  increases  rapidly  for  the  smallest  diameter  particles.  The  derivative  is 
negative  throughout  the  diameter  space  indicating  a  decrease  in  the  flux  as  diameter  increases. 
Cases  I  and  n  (Pig.  22)  show  the  similar  characteristics. 
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Figure  22.  Flux  model  sensitivity  with  respect  to  debris  diameter  (Cases  I  &  H). 


The  flux  is  extremely  sensitive  to  changes  in  diameter  for  the  smaller  particles,  especially  beginning 
where  l3F/9dl  >  1  at  d  <  0.05  cm.  This  is  due  to  the  nature  of  the  curve  fit,  Fj,  through  the  data 
(Fig,  1).  The  function  F(d,  h,  i,  t,  S)  predicts  an  infinite  flux  as  diameter  approaches  zero 
(discussed  in  Subsection  3.4).  Because  of  this  behavior,  it  is  recommended  that  the  model  be 
limited  to  the  range  of  diameters,  10-  <  d  <  lO^  cm,  where  the  basis  is  measured  data. 

6.2.2  3F/9h  versus  h 

The  variation  of  3F/3h  versus  altitude  is  shown  in  Figure  24  for  the  two  base  line  cases.  The 
magnitude  of  aF/ah  is  very  small,  i.e.,  SF/ahl «  1;  thus,  while  there  is  change,  the  predicted  flux 
is  relatively  insensitive  to  changes  in  altitude.  From  Figure  24,  it  is  seen  that  the  sensitivity  of 
aF/ah  increases  through  the  lower  altitudes  and  peaks  at  -500  km  where  much  of  the  data  base 
lies.  Then  the  sensitivity  decreases  with  increasing  altitude. 
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(Impacts 


Figure  24,  Flux  model  sensitivity  with  respect  to  orbital  altitude  (Cases  I  &  H). 


The  predicted  flux  is  relatively  insensitive  to  changes  in  altitude,  especially  as  altitude  increases. 
The  lack  of  flux  changes  with  altitude  beyond  1000  km  suggests  that  the  model  be  limited  to 
altitudes  <1000  km.  The  sensitivity  increases  through  the  lower  altitudes,  up  to  -500  km,  then 
decreases  (becoming  less  sensitive)  as  altitude  increases  beyond  500  km. 

6.2.3  BF/di  versus  i 

The  variation  of  3F/ai  versus  orbital  inclination  is  shown  in  Figure  25  for  the  two  base  line  cases. 
The  curves  shown  reflect  the  nature  of  'F.  The  expression  3F/3i  is  equal  to  zero  at  80, 90,  and 
100  deg  where  'P  has  local  maximums  and  a  minimum.  The  function  is  relatively  insensitive  to 
changes  in  orbital  inclination,  i.e.,  3F/3i «  1. 

The  flux  prediction  is  relatively  insensitive  to  changes  in  orbital  inclination  and  mirrors  the 
variation  in  T'  with  respect  to  i. 


45 


1 


Figure  25.  Flux  model  sensitivity  with  respect  to  orbital  inclination  for  Cases  I  &  n, 
where  d  =  1  cm,  h  =  500  km,  p  =  0.05,  and  q  =  0.02. 

6.2.4  3F/3t  versus  t 

The  variation  of  dF/dt  versus  time  (yrs)  is  shown  in  Figure  26  for  the  two  base  line  cases.  The 
curves  show  an  increase  in  sensitivity  with  time,  but  the  magnitude  (10-  '^)  is  extremely  small,  i.e., 

aF/at «  1. 

The  flux  prediction  sensitivity  increases  with  time,  but  the  overall  magnitude  is  extremely  small. 
This  increased  sensitivity  is  negligible  over  the  period  of  time  for  which  the  model  is  being 
considered  valid. 

6.2.5  aF/aS  versus  S 

The  variation  of  dF/dS  versus  the  measure  of  the  solar  activity,  S,  is  shown  in  Figure  27  for  the 
two  base  line  cases.  The  magnimde  of  the  sensitivity  is  extremely  small,  i.e.,  dF/dS  «  1.  As 
seen  in  Figure  27,  laF/aSI  reaches  a  maximum  at  approximately  S  =  140,  which  is  the  value  used 
to  normali2e  S  in  ^>i. 
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The  flux  prediction  sensitivity  with  respect  to  S  is  extremely  small  and  on  the  same  order  of 
magnitude  as  the  previous  sensitivities. 


6.2.6  9F/3p  and8F/3q 


The  variations  of  3F/9p  versiB  p  and  3F/3q  versus  q  are  shown  in  Figure  28  for  Case  n  only 
(t  =  1995).  Case  I  is  trivial  since  both  gi  and  g2  equal  one  for  t  =  1988,  i.e.,  no  change  in  the 
sensitivities  with  either  p  or  q.  As  can  be  seen  in  Figure  28, 3F/9p  is  constant  with  respect  to 
changes  in  p.  This  is  due  to  the  linear  nature  of  g2.  The  function  3F/Bq,  however,  changes  with  q 
and  increases  when  q  is  increased,  reflecting  the  compounded  nature  of  g2.  This  points  out 
another  difference  between  compounded  growth  (gi)  or  linear  growth  (g2),  discussed  in 
Subsection  4.3.  The  flux  prediction  sensitivity  does  not  change  for  linear  growth,  while  its 
sensitivity  increases  with  compounded  growth  (increases  with  increasing  q).  The  relative  order  of 
magnitude  of  9F/9q  10-  ranks  this  sensitivity  as  second  to  3F/3d. 

The  flux  prediction  sensitivity  does  not  change  with  the  large-particle  growth  rate,  p,  and  is 
increasingly  sensitive  to  increases  in  the  small-particle  growth  rate,  q. 


Growth  Rates,  p  and  q 

Figure  28.  Flux  model  sensitivity  with  respect  to  particle  growth  rates  (Case  H). 


6.2.7  Combined  Sensitivities  for  3^  F/3h  3d.  3^  F/^h  and  ^2  F/;)d  ^g. 

Combined  or  composite  sensitivities  are  shown  in  Figures  29  to  31  for  Case  H.  Because  of  the 
nature  of  Equation  6,  the  sensitivities  presented  are  given  by  the  following  relationships: 
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Figure  29  shows  the  flux  model  sensitivity  versus  altitude  and  diameter.  As  seen  in  this  figure,  the 
flux  model  is  most  sensitive  for  the  smaller  diameters  at  the  lower  altitudes  and  becomes  extremely 
insensitive  at  the  higher  altitudes  for  the  larger  particles. 


The  flux  model  sensitivity  versus  altitude  and  inclination  given  in  Figure  30  shows  the  flux  is 
sensitive  to  variation  in  orbital  inclination  for  only  the  lower  altitudes,  h  <  1000  km.  Above 
1000  km,  the  model  is  insensitive  to  any  changes  in  altitude  or  orbital  inclination  and  a  uniform 
debris  cloud  is  predicted  above  1000  km. 


The  variation  of  the  flux  model  sensitivity  with  respect  to  diameter  and  small-particle  growth  rate 
is  shown  in  Figure  31.  Clearly  the  dominant  sensitivity  is  with  respect  to  diameter.  The  results  of 
these  sensitivities  are  consistent  with  the  results  previously  presented.  The  flux  model  is  most 
sensitive  to  changes  in  diameter,  especially  as  diameter  decreases.  The  flux  model  is  least  sensitive 
to  changes  in  altitude,  predicting  a  constant  debris  flux  beyond  1000  km. 
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Figure  29.  Flux  model  sensitivity  with  respect  to  altitude  and  diameter  (Case  II), 
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Figure  31.  Flux  model  sensitivity  with  respect  to  diameter  and 
small-particle  growth  rate  (Case  II). 


6.3  SENSrnVITY  ASSESSMENT 

The  sensitivities  of  the  flux  model  have  been  analyzed  with  respect  to  the  independent  variables, 
debris  diameter  (d),  orbital  altitude  (h),  orbital  inclination  (i),  time  (t),  solar  activity  (S),  and  the 
large-  and  small-particle  growth  rates  (p  and  q).  The  summary  results,  depicted  in  Figure  32, 
show  the  relative  magmtudes  of  these  sensitivities.  One  can  see  the  flux  prediction  is  most 
sensitive  to  changes  in  debris  diameter,  for  the  small  diameter  debris.  On  the  basis  of  the  singular 
behavior  of  F  as  diameter  decreases  and  the  range  of  diameters  from  the  data  base,  it  is 
recommended  the  model  be  restricted  to  10-  ^  <  d  <  KP  cm. 


Next  in  importance  is  the  sensitivity  of  F  to  the  assumed  growth  rates  for  small  particles,  q.  The 
sensitivity  increases  as  the  growth  rate  increases.  The  sensitivities  with  respect  to  the  large-particle 
growth  rate,  p,  and  inclination  are  next  in  relative  importance.  The  sensitivity  with  respect  to  p  is 
constant  The  sensitivity  with  respect  to  i  at  500-km  mirrors  Y .  (Note  that  for  comparison  the 
orbital  inclination  sensitivity  presented  in  Figure  32  has  the  units  of  flux/radian.)  The  sensitivities 
with  respect  to  time  and  solar  activity  are  next  in  importance  but  are  much  less  in  magnitude  than 
the  previous  variables.  The  least  sensitive  changes  by  an  order  of  magnitude  are  with  respect  to 


f, 


0  0.2  0.4  0.6  0.8  1 

Normalized  Variable  (d,h,i,t,S,p,q) 


Figure  32.  Flux  model  sensitivities,  summary  comparisons  (Case  II).  (Baseline  values: 
h  =  500  km,  d  =  1  cm,  i  =  47  deg,  t  =  1995,  S  =  90,  p  =  5%,  q  =  2%.) 

altitude.  This  suggests  that  the  model  does  not  account  for  variations  with  altitude  in  a  realistic 
manner.  In  fact,  F  becomes  nearly  constant  as  altitude  increases  beyond  nominally  800  to 
1000  km  (3F/9h  0).  Recall  that  Equation  6  does  not  reproduce  the  published  data  as  was 

shown  in  Figure  3  for  these  altitudes.  This  insensitivity  and  weakness  in  the  model  prediction  with 
respect  to  altitude  emphasize  its  limitations  and  suggest  the  model  be  used  at  altitudes  <1000  km 
rather  than  what  is  typically  thought  of  as  LEO  (h  <  2000  km).  In  summary,  the  rank  order  of  the 
flux  model  sensitivities  is  as  follows: 

Parameter  Relative  Order  of  Magnitude/Comments 

•  d  (diameter)  ~10  ^  at  d  =  0.001  cm 

~1  at  d  =  0.05  cm 
~10-9atd=  100  cm 

Flux  prediction  becomes  undefined  as  d  — »  0. 

•  q  (smaU-partrcle  growth  rate)  ~10‘  slightly  increasing  with  increasing  q 
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•  i  (orbital  inclination)  and 

p  (large-particle  growth  rate) 

•  t  (time)  and  S  (solar) 

•  h  (altitude) 


~10'  5  to  10'  sensitivity  with  respect  to  i  reflects 
'P(i)  variation;  sensitivity  with  respect  to  p  is  constant 

~10-7tol0-8 


~10*  8  (h  =  500  km)  to  10'  (h  =  2000  km),  predicts 
constant  flux  for  h  >  10(X)  km 
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7.0  UNCERTAINTY  ANALYSIS 


This  section  presents  an  uncertainty  analysis  of  the  NASA  90  model.  The  uncertainty  analysis  is  a 
study  of  how  the  root-mean-square  (RMS)  error  in  the  data  or  curve  fit  affects  the  results  obtained 
from  the  use  of  the  model.  The  uncertainty  analysis  approach  is  described  in  Subsection  7.1  while 
Subsection  7.2  shows  the  results.  Subsection  7.3  presents  the  summary  assessment. 

7.1  PROPAGATION  OF  ERRORS 

Most  empirically  based  models  will  have  some  inherent  error  in  the  results  because  of  the 
uncertainties  in  the  measured  values  and  because  of  the  RMS  difference  of  the  curve  fit  to  the 
These  differences  can  be  due  to  either  the  natural  variation  of  parameters  measured  in  space  and 
time  or  in  the  model  (curve  fit)  representation  of  the  data.  The  values  of  the  differences,  or  errors, 
can  be  estimated  with  some  level  of  confidence  from  the  known  or  quantifiable  characteristics  of 
the  data  or  the  curve  fit  For  example, 

±  ej  =  estimated  uncertainty  in  the  parameter  j  at  some  level  of  confidence  (e.g.,  90% 
confidence).  Note  that  the  +/-  values  may  not  be  necessarily  identical. 

For  the  NASA  90  model,  the  overall  error,  ±  ep,  in  F(d,  h,  i,  t,  S)  is  due  to  the  uncertainties  in  the 
variables  used  to  determine  F.  If  one  assumes  that  each  of  the  variables  are  uncorrelated,  the 
contribution  of  the  individual  uncertainties  to  ep  may  be  estimated  by  the  use  of  a  propagation  of 
error  analysis  as  follows: 


where 

Ep  =  Uncertainty  in  F,  e.g.,  F  ±  Ep 

Ed  =  Uncertainty  in  flux  with  respect  to  particle  diameters 

Eh  =  Uncertainty  in  flux  with  respect  to  altitude 

Ej  =  Uncertainty  in  flux  with  respect  to  inclination 

Et  =  Uncertainty  in  flux  with  respect  to  time 

Ep  =  Uncertainty  in  flux  with  respect  to  large-particle  growth 
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£q  =  Uncertainty  in  flux  with  respect  to  small-particle  growth 
Es  =  Uncertainty  in  the  level  of  solar  activity 


Each  of  the  individual  uncertainties  is  estimated  at  the  same  level  of  confidence.  In  this  study 
NASA's  estimates  of  the  individual  uncertainties  as  presented  in  Subsection  4.6  were  used.  This 
led  to  the  modification  of  Equation  27  as  follows: 


Ef  =  ± 


(efm)^ +(eFH)^  + 


dF  f 
—  es  + 
\as  / 


1/2 


(28) 


where 


epM  =  Uncertainty  in  F  because  of  statistical  and  measurement 
uncertainties  primarily  as  a  function  of  debris  diameter 

Efh  =  Uncertainty  in  F  with  respect  to  altitude  due  in  part  to  the  difficulty 

in  determining  F  for  debris  in  highly  elliptical  orbits 

Equation  28  represents  the  parameters  that  were  felt  to  affect  significantly  the  uncertainty  in  the  flux 
prediction,  Equation  6.  Note  that  aF/ap,  dF/dq,  and  aF/aS  are  given  by  Equations  21  through  23. 
The  estimated  uncertainties  used  for  this  study  are  shown  in  Table  5. 

Rationale  for  Choices:  Equation  28  was  used  to  analyze  the  results  of  the  error  propagation 
because  it  was  convenient  to  use  the  uncertainties  reported  by  NASA  (Ref.  2).  The  terms  shown 
in  Equation  27  are  equivalent  to  those  in  Equation  28.  The  uncertainty  in  flux  measurements,  Efm 
and  the  uncertainty  with  respect  to  the  altitude  distribution,  Epu,  in  Equation  28  represent  the 
measurement  uncertainties  given  as  the  first  four  terms  of  Equation  27.  Because  Efm  and  Efh  are 
reported  directly  by  NASA  (Ref.  2),  there  is  no  need  to  attempt  to  estimate  Ed ,  Eh ,  Ei ,  and  Et 
directly.  The  remaining  uncertainties  in  knowledge  of  the  actual  particle  growth  and  solar  activity, 
Ep,  Eq,  and  Eg,  are  estimated  by  NASA  (Ref.  2)  and  were  included  because  the  state  of  particle 
growth  has  not  yet  been  verified  and  because  predicting  future  maximum  levels  of  solar  activity 
remains  uncertain  (Ref.  17). 

An  independent  assessment  of  the  NASA  estimated  uncertainties  was  not  conducted  as  part  of  this 
task.  Such  an  assessment  may  be  performed  by  establishing  uncertainties  in  each  variable  in 
Equation  27. 
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Table  5.  Propagation  of  error  analysis  uncertainties. 


Parameter 

Estimated  Uncertainty 

Flux  Uncertainty 

(90%  Confidence) 

•  Flux  Measurements: 

d  >  10  cm 

iepM  ='0.5F 

0.05  <  d  <  10  cm 

+eFM  =  2F 
-epM  =  (2/3)  F 

d  <  0.05  cm 

+£pM  =  F 
-epM  =  0.5  F 

•  Flux  Measurements:  Altitude 

Distribution 

d  >  10  cm 

+eFH  =  F 

-EpH  =  0.5  F 

1  <  d  <  10  cm 

i£pH  =  (same  as  ±£fm  for  0.05  <  d< 

10  cm) 

d<  1  cm 

lepH  =  4F  for 

100  <  h  <  300  km 

±eFH  =  0  for  300  <  h  <  700  km 

±£FH  =4F, 
for  700  ^  h  <  900  km 

±Efh  =  9F  for  900  <  h  <  1100  km 

• 

• 

±5F  for  each  200-km  increment 
in  altitude 

• 

• 

• 

±epH  =  34F  for  1900  <  h  <  2000  km 

•  Large-Particle  Growth 
(p  =  0.05) 

+£  p  =  0.05 

-£  p  =  0.01 

•  Small-Particle  Growth 
(q  =  0.02) 

+£  q  =  0.18 

-£  q  =  0.02 

•  Solar  Activity 

+£  s  =  250  -  S 
-£  s  =  S  -  70 
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7.2  UNCERTAINTY  ANALYSIS  RESTTTTS 


This  subsection  discusses  the  uncertainties  given  in  Equation  28.  The  Em  and  Efh  uncertainties 
will  be  analyzed  in  Subsections  7.2.1  and  7.2.2.  The  contribution  to  +eF ,  +eF  *>  of  two  of  the 
uncertainties,  +£fm  and  +£fh»  are  plotted  in  Figure  33  as  a  function  of  d  and  h.  The  flux, 
F(h,d,i,t,p,q),  has  been  factored  out,  i.e.,  normalized  error,  +£f  *■ 

0  0  1/2 

+  £p  oc  (cjF  +  C2F^)  (29) 

or 


+  ep*  =  F  (ci+C2)^^^ 

^  F 


(30) 


where  Cj  and  C2  are  the  multipliers  of  F  (Table  5)  that  are  used  to  calculate  +£fm  and  +£fh  for 
varying  diameters  and  altitudes. 

Figure  33  shows  that  the  normalized  error,  h-Ef  *,  is  greater  than  a  factor  of  5  for  altitudes 
>700  km  and  particle  diameters  <1  cm.  For  altitudes  >1000  km  and  particle  diameters  <1  cm,  +£f 

*  is  greater  than  one  order  of  magnitude.  These  observations  will  help  in  visualizing  the  overall, 
uimormalized  error  shown  in  Figures  34-39. 


Figure  33.  Flux  model  imcertainty,  -h£fm  and  +eFH  components,  versus  diameter  and 
altitude. 
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The  graphical  results  of  the  absolute  uncertainty,  ±£f,  are  presented  in  Figures  34-36  and 
Figures  37-39  for  the  AF/DoD  orbits  of  interest  identified  in  Section  5.0.  Figures  34-36  present 
the  flux  versus  the  diameter  for  orbital  altimdes  of  500, 1000,  and  1500  km  and  orbital  inclinations 
of  30, 70,  and  95  deg,  respectively.  The  debris  diameter  is  varied  between  lO-'*  to  10^  cm. 

Figures  36-38  present  the  flux  versus  altitude  for  diameters  of  0.05, 0.l,  1.0,  and  10  cm  and 
orbital  inclinations  of  30, 70,  and  95  deg.  The  orbital  altitude  is  varied  between  100  to  2000  km. 
In  all  figures,  the  upper  curve  represents  F  +  ep,  the  middle  curve  represents  the  nominal  value  of 
F  and  the  lower  curve  represents  F  -  ep.  The  uncertainties,  +ep  and  -ep,  were  determined  using 
Equation  28  and  Table  5. 

7.2.1  Flux  versus  Diameter 

The  flux  versus  diameter  results  are  presented  in  Figures  34-36  (a,  b,  c)  for  1995  (S  =  90). 

Figure  34  presents  the  results  for  i  =  30  deg.  Figure  35  for  i  =  70  deg,  and  Figure  36  for 
i  =  95  deg  for  orbital  altimdes  of  (a)  500  km,  (b)  1000  km,  and  (c)  1500  km.  An  examination  of 
these  figures  shows  that  the  uncertainty  increases  with  increasing  altimde  and  decreases  with 
increasing  diameter.  The  best  results  (those  with  the  least  uncertainty)  are  for  an  orbital  altimde  of 
500  km  and  diameters  >10  cm.  This  is  expected  because  much  of  the  impact  data  base  comes  from 
orbital  altimdes  of  500  km  or  less  (Solar  Max).  The  best  measurements  with  respect  to  size  come 
from  USSPACECOM  data,  which  have  an  estimated  measurement  resolution  of  diameters  down  to 
10  cm. 


Typical  of  the  general  trends,  the  flux  uncertainty  varies  by  a  factor  of  2  to  4  for  the  lower  altimdes 
and  by  an  order  of  rnagnimde  for  the  higher  altitudes.  The  major  contributors  to  this  overall 
uncertainty  can  be  identified  by  examining  the  individual  terms  of  Equation  28.  For  example,  at 
i  =  30  deg,  h  =  1000  km,  d  =  0.1  cm  (Fig.  33),  and  F  =  0.00407  impacts/m^  •  yr. 


ep  =  ± 


(epm)^  +  (eph)^  + 


1/2 


=  [(8.15  X  10-3)2  +  (3.67  X  10-2)2  +  (2.21  x  10-2)2+  (5.03  x  10-3)2+  (1.49  x  10-5)2]l/2 
=  3.79  (10-2)  thus  [(F  +  £f)/F]  =  10.3 
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The  flux  prediction  is  dominated  by  the  first  two  terms  which  represent  the  uncertainties  in  the  flux 
measurement  capabilities.  This  emphasizes  earlier  conclusions  that  the  accuracy  of  this  model  can 
be  increased  significantly  through  a  reduction  in  the  measurement  uncertainties.  The  next  term  to 
dominate  is  the  uncertainty  in  the  small-particle  growth  (recall  that  this  example  is  for 
d  =  0.1cm). 


(b)  h  =  1000  km,  t  =  1995,  and  S  =  90. 

Figure  34.  Flux  model  uncertainty  versus  diameter  for  i  =  30  deg. 


59 


Flux  (Impacts/m  ^yr) 


I  I  I  I  I  llllj  T  I  I  I  Mll|  I  III  I  j  I  I  I  I  llllj  I  I  I  I  II  II  j  I  I  I  I  II  II  [  I  I  I  I  II  ll| 

0.0001  0.001  0.01  0.1  1  10  100  1000 

Diameter  (cm) 

(c)  h  =  1500  km,  t  =  1995,  and  S  =  90. 

Figure  34.  Concluded. 
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(a)  h  =  500  km,  t  =  1995,  and  S  =  90. 

Figure  35.  Rux  model  uncertainty  versus  diameter  for  i  =  70  deg. 
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(c)  h  =  1500  km,  t  =  1995,  and  S  =  90. 
Figure  35.  Concluded. 
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(c)  h  =  1500  km,  t  =  1995,  and  S  =  90. 

Figure  36.  Concluded. 

7.2.2  Flux  Versus  Altitude 

The  flux  versus  altitude  results  are  presented  in  Figures  37-39  (a,b,c,d)  for  1995  (S  =  90). 

Figure  37  presents  the  results  for  i  =  30  deg.  Figure  38  for  i  =  70  deg,  and  Figure  39  for 
i  =  95  deg  for  orbital-debris  particle  diameters  of  (a)  0.05  cm,  (b)  0.1  cm,  (c)  1.0  cm,  and  (d) 

10  cm.  Examination  of  these  figures  shows  the  least  uncertainty  for  altitudes  around  500  km, 
with  the  uncertainty  increasing  as  one  moves  away  from  this  altitude. 

The  flux  uncertainty  can  vary  by  a  factor  of  2  for  the  large  diameters  and  by  an  order  of  magnitude 
for  the  smaller  diameters.  The  uncertainty  in  F,  F  +  ep,  increases  as  the  altitude  approaches 
2000  km.  The  major  contributors  to  these  levels  of  uncertainty  can  again  be  identified  by 
examining  the  terms  of  Equation  28.  For  example  (Fig.  37a),  at  i  =  30  deg,  h  =  1000  km, 
d  =  0.05  cm,  and  F  =  0.023  impacts/m^  •  yr. 


=  ±  ( 


eF  =  ±  (eFM)2+(epH)2  + 


Up 


Ep 


.2  / 

+ 


Uq  /  Us 


es 


2-11/2 


=  [(2.3  X  10-2)2  +  (2.07  x  10-1)2  +  (2.21  x  10-7)2  +  (2.85  x  10-2)2  +  (g  41^  10*5)2]  1/2 
=  2.11  (10-1)  aus.  [(F  +  eF)/F]  =  10.1 
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As  in  the  case  for  flux  versus  diameter,  the  dominant  terms  are  the  uncertainties  in  the  flux 
measurements  and  the  small-particle  growth.  The  uncertainties  with  respect  to  large-particle 
growth  and  the  level  of  solar  activity  are  negligible  because  the  choice  of  diameters  for  this  example 
(d  =  0.05  cm,  small  particle)  and  choice  of  altitude  where  solar  effects  are  minimal. 


(a)  d  =  0.05  cm,  t  =  1995,  and  S  =  90. 


(b)  d  =  0.1  cm,  t  =  1995,  and  S  =  90. 

Figure  37.  Flux  model  uncertainty  versus  altitude  for  i  =  30  deg. 
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(b)  d  =  0.1  cm,  t  =  1995,  and  S  =  90. 

Figure  38.  Flux  model  uncertainty  versus  altitude  for  i  =  70  deg. 
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Flux  ampacts/m2*yr)  Flux  (Impacts/m  .yr) 


(d)  d=  10cm,t=  1995,8  =  90. 
Figure  38.  Concluded. 
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»,irapacis/m“*  yr) 


Figure  39.  Concluded 
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7.3  UNCERTAINTY  ASSESSMENT 


Table  6  summarizes  the  imcertainties  analyses  presented  for  an  orbital  inclination  of  95  deg  and 
t  =  1995  (S  =  90)  for  the  altitudes  and  diameters  shown.  Figure  40  presents  the  overall  model 
uncertainty  versus  diameter  and  altitude  for  1995  and  an  orbital  inclination  of  30  deg. 

From  Table  6,  except  for  an  altitude  of  500  ±  200  km  and  diameters  ^10  cm,  the  model 
uncertainties  are  at  least  plus  or  minus  an  order  of  magnitude  and  increase  as  altitude  approaches 
2000  km.  For  aU  altitudes  and  debris  particle  sizes,  the  overall  uncertainty,  ep,  is  dominated  by  the 
uncertainty  in  the  measurements,  Efm  and  epH-  The  smaU-particle  growth  uncertainty,  (3F/3q)eq,  is 
significant  for  the  small  particles,  d  <  10  cm,  while  the  large-particle  growth  uncertainty, 
(3F/3p)ep,  gains  importance  as  particle  size  increases  beyond  d  ^  10  cm.  The  relative  influence  of 
solar  activity,  as  measured  by  the  value  of  S,  increases  as  altitude  decreases  (for  a  given  size 
particle)  for  altitudes  <1000  km.  Above  these  altitudes  the  influence  of  the  atmosphere  is 
negligible,  regardless  of  the  value  of  S,  on  the  overall  flux  uncertainty. 

Table  6.  Orbital-debris  flux  uncertainty 
(i  =  95deg,t=  1995,  S=90). 


h(km)  d(cm)  epM 


epH 


ep  F-t-gp 
F 


500 

0.01 

1.51 

7.55(10-1) 

♦ 

1.86 

1.21 

2.79 

2.85 

0.10 

9.55(10-3) 

2.39(10-3) 

* 

5.90(10-3) 

3.84(10-3) 

1.21(10-2) 

3.53 

1.00 

4.16(10-5) 

1.04(10-5) 

.  * 

2.41(10-5) 

1.67(10-5) 

5.19(10-5) 

3.50 

10.00 

1.34(10-5) 

2.69(10-5) 

6.62(10-7) 

1.64(10-7) 

2.16(10-5) 

3.76(10-5) 

2.40 

100.00 

2.31(10-7) 

4.63(10-7) 

1.20(10-7) 

* 

3.72(10-7) 

6.48(10-7) 

2.40 

1000 

0.01 

2.17 

19.5 

* 

2.68 

19.8 

10.10 

0.10 

1.37(10-2) 

6.18(10-2) 

* 

8.48(10-3) 

* 

6.38(10-2) 

10.30 

1.00 

5.98(10-5) 

2.69(10-‘l) 

* 

3.46(10-5) 

* 

2.78(10-3) 

10.30 

10.00 

1.93(10-5) 

3.86(10-5) 

9.52(10-7) 

2.53(10-7) 

DC 

4.43(10-5) 

2.15 

100.00 

3.33(10-7) 

6.65(10-7) 

1.72(10-7) 

* 

♦ 

7.63(10-7) 

2.15 

1500 

0.01 

2.17 

52.1 

* 

2.68 

* 

52.2 

25.10 

0.10 

1.37(10-2) 

1.65(10-1) 

* 

8.49(10-3) 

* 

1.66(10-1) 

25.10 

1.00 

5.99(10-5) 

7.18(10-4) 

* 

3.47(10-5) 

* 

7.22(10-4) 

25.10 

10.00 

1.93(10-5) 

3.87(10-5) 

9.53(10-7) 

2.36(10-5) 

* 

4.43(10-5) 

2.15 

100.00 

3.33(10-7) 

6.66(10-7) 

1.73(10-7) 

* 

* 

7.64(10-7) 

2.15 

Note:  *  Small  Relative  to  Other  Terms 
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Figure  40  shows  the  regions  where  the  flux  model  gives  the  best  prediction.  The  uncertainty 
factor,  [(F  +  eF)/F],  versus  diameter  and  altitude  for  i  =  30  deg,  t  =  1995,  and  S  =  90.  The  most 
precise  predictions  are  at  an  altitude  of  500  km  (±200  km)  or  for  the  large  particles,  d  >  5  cm. 
These  results  are  consistent  with  the  greatest  availability  of  data.  Solar  Max  and  for  these  size 
particles,  USSPACECOM.  The  uncertainty  increases  as  particle  size  decreases  or  as  altitude 
increases  away  from  these  nominal  values.  The  overall  uncertainty  for  these  latter  regions  is 
shown  to  be  at  least  plus  or  minus  an  order  of  magnitude. 


As  shown  by  the  examples  in  Subsection  7.2  and  as  depicted  in  Table  6  and  Figure  40,  the  overall 
uncertainty  of  this  model  is  consistently  dominated  by  the  uncertainties  in  the  measured 
environment,  efm  efh-  If  one  could  improve  the  space  debris  measurements,  i.e.,  reduce 
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these  uncertainties,  especially  epH  with  respect  to  altitude,  then  the  overall  model  uncertainty  could 
be  reduced.  This  assumes  that  the  model  accurately  predicts  the  environment  to  begin  with  and  that 
the  improved  measurements  are  consistent  with  the  model  predictions  (there  exists  the  possibility 
that  improved  measurements  could  invalidate  the  model).  With  improved  measurements,  reduction 
in  the  uncertainties  are  possible  for  both  the  large-  and  small-particle  growth,  Ep  and  Eq.  For 
example,  an  order  of  magnitude  reduction  in  ejm  and  efh  would  result  in  the  following  for  a 
particle  size  of  0.1  cm  at  30-deg  orbital  inclination  and  altitudes  of  1000  and  1500  km  (Table  7). 


Table  7.  Improved  flux  predictions. 


dfcml 

hflanl 

£fm 

£fh 

£f 

F 

(F+ef)/F 

Current 

0.1 

1000 

8.15(10-3) 

3.67(10-2) 

3.79(10-2) 

o 

-j 

o 

1 

10.3 

1500 

8.16(10-3) 

9.79(10-2) 

9.84(10-2) 

4.08(10-3) 

25.1 

Improved 

0.1 

1000 

8.15(10-4) 

3.67(10-3) 

6.28(10-3) 

4.07(10-3) 

2.54 

1500 

8.16(10-4) 

9.79(10-3) 

1.10(10-2) 

4.08(10-3) 

3.71 

Improving  the  measurement  uncertainties  by  an  order  of  magnitude  can  reduce  the  overall  model 
uncertainty  to  a  consistent  factor  for  all  altitudes  and  sizes  and  to  values  less  than  an  order  of 
magnitude. 
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8.0  CONCLUSION 


The  NASA  90  orbital-debris  environment  prediction  model  has  been  assessed.  The  model 
represents  a  smooth  curve  fit  to  the  data.  These  data  consist  of  USSPACECOM  radar  data, 
MIT/LL  ETS  and  GEODSS  optical  data,  returned  spacecraft  surface  impact  data  (Solar  Max),  and 
recent  radar  data  from  Aricebo  and  Goldstone.  Special  purpose  analytic/computational  models 
were  used  by  NASA  to  validate  future  traffic  and  growth  studies.  The  USSPACECOM  data  were 
used  primarily  to  establish  altitude,  velocity,  and  solar  and  orbital  inclination  effects  for  debris 
particles  10  cm  or  greater  in  diameter.  Particle-size  distribution  was  established  by  examining  all 
data  sets,  but  primarily  relied  on  impact  studies  of  returned  spacecraft  surfaces  for  the  smaller- 
debris  particles. 

8.1  MODEL  DEVELOPMFNT 

NASA  believes  that  the  NASA  90  model  is  valid  only  for  spacecraft  in  LEO  and  through  the  year 
2010.  The  validity  of  the  model  rests  on  the  fundamental  assumptions  made  in  the  derivation  of 
this  model.  Other  than  the  assumptions  concerning  the  vaHdity  of  the  data  or  the  specific  uses  of 
those  data  in  the  model  formulation,  perhaps  the  most  fundamental  or  key  observation  concerning 
this  model  is  with  respect  to  that  data  base.  The  model  was  derived  to  fit  the  data  that  represent  a 
historic  record  of  man's  first  space  operations.  If  these  operations  should  change  dramatically, 
then  the  basis  for  this  model  may  no  longer  be  valid.  For  example,  the  model  assumes  certain 
growth  rates,  breakup  rates,  and  has  inherent  debris  concentrations,  ¥(1),  built  into  it  If  any  of 
these  drastically  change,  the  model  may  not  be  adjustable  to  account  for  these  changes  (nor  may  the 
model  be  used  to  assess  potential  space  architectures  where  these  ground  rules  change 
significantly).  Even  now  there  are  significant  altitudes  where  the  smoothed  NASA  model  does  not 
fit  the  debris  environment  by  one  to  two  orders  of  magnitude  because  of  breakups  at  these  altitudes 
and  inclinations.  The  model  accuracy  is  fundamentally  tied  to  the  accuracy  of  the  data  base.  The 
data  base  is  most  accurate  for  altitudes  on  the  order  of  500  km  and  debris  particles  >10  cm  in 
diameter.  The  model  accuracy  can  only  be  improved  by  an  increase  in  the  accuracy  of  the 
assuming  any  new  or  improved  data  can  be  incorporated  easily  into  the  model.  For  example,  the 
function  H(d)  was  developed  to  account  for  improved  knowledge  of  the  environment  for  particles 
>8  cm  (although  extrapolated  to  smaller  sizes).  Future  fixes  may  not  be  as  easily  developed  or 
incorporated  to  improve  the  model's  performance.  Perhaps  the  most  controversial  and  uncertain 
part  of  the  model  is  in  predicting  the  future.  Tliere  simply  is  not  enough  of  a  historic  record  to  see 
the  future  trends  or  growth  rates  to  accurately  project  whether  these  are  indeed  linear  or 
compounded  growth  curves. 


73 


The  most  unsatisfactory  parts  of  this  model  are  the  collision-velocity  impact  distributions,  velocity 
direction,  and  the  particle-density  estimates.  The  collision- velocity  impact  distribution  is  a 
complicated  function  of  orbital  inclination  and  has  no  discernible  link  to  physical  reality  in  its 
presented  form.  The  direction  of  impact  model  assumes  circular  orbits  only  and  does  not  allow  for 
elliptical  orbit  out-of-plane  impacts.  This  may  not  be  a  bad  assumption  as  a  first  approximation 
since  most  LEO-region  debris  orbits  will  be  nearly  circular.  The  particle-density  function  assumes 
spherical  particles  and  has  little  data  to  support  its  formulation,  especially  for  the  smaller  sizes. 
Each  of  these  expressions  are  restricted  due  to  limited  data.  Improving  these  expressions,  or 
knowledge  of  these  environments,  will  only  come  with  improved  measurements  and  experiments. 
The  user  is  cautioned  that  these  particular  expressions  in  the  NASA  90  model  are  particularly 
weak. 

8.2  SENSmVTTY  ANALYSTS 

The  NASA  90  orbital-debris  flux  prediction  model, 

F  (d,  h,  i,  t,  S)  =  H  (d)  ^  (h,  S)  ^  (i)  [Fi  (d)  gi  (t)  +  F2  (d)  g2  (t)] 

has  been  shown  to  be  most  sensitive  to  changes  in  the  orbital-debris  particle  diameter,  d.  The 
changes  in  flux  are  most  sensitive  to  diameters  <0.05  cm,  with  the  overall  sensitivity  increasing  as 
particle  size  decreases.  This  sensitivity  reflects  the  nature  of  the  curve  fit  predictions,  F^  (d)  and 
F2  (d).  The  flux  model  is  relatively  insensitive  to  the  large  diameters  >1  cm. 

Next  in  relative  importance  is  the  model  sensitivity  with  respect  to  the  small-particle  growth  rate,  q, 
followed  by  the  sensitivities  with  respect  to  large-particle  growth,  p,  and  inclination,  i.  The  model 
sensitivity  increases  as  q  increases,  which  reflects  the  nature  of  the  assumed  compounded  growth, 
gi(t),  for  the  small  particles.  Changes  in  the  value  of  p  do  not  influence  the  model  sensitivity, 
reflecting  the  nature  of  the  assumed  linear  growth,  g2(t),  for  the  large  particles.  This  difference  in 
sensitivity  illustrates  the  distinction  between  whether  debris  growth  is  assumed  to  be  linear  or 
compounded.  Neither  linear  nor  compounded  growth  for  the  large-particles  is  substantiated  by  the 
data.  There  is  essentially  no  data  to  base  the  small-particle  growth  rate.  If  the  actual  growth  rates 
are  found  to  be  compounded  for  p  and  q,  then  the  choice  of  the  growth  rate  (5  percent  versus 
10  percent)  will  influence  the  resulting  prediction.  If  the  growth  is  linear,  then  the  actual  value  of 
the  growth  rate  becomes  less  important.  The  sensitivity  with  respect  to  orbital  inclination  reflects 
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the  'P(i)  distribution  defining  regions  of  increased  and  decreased  sensitivity  as  the  inclination  is 
varied  for  altitudes  <1000  km. 

The  flux  model  sensitivity  to  time  (t)  and  level  of  solar  activity  (S)  are  relatively  equal  in 
importance.  The  least  sensitive  parameter  is  altitude  (h).  Beyond  800  to  1000-km  altitude,  the  flux 
model  predicts  a  cloud  of  orbital  debris  that  remains  constant  with  increasing  altitude.  This  does 
not  fit  physical  intuition  nor  does  it  comply  with  the  data  as  shown  in  Figure  3.  In  Figure  3,  the 
debris  flux  environment  peaks  at  800, 1000,  and  1500  km  because  of  known  breakups  at  these 
altitudes,  which  the  model  clearly  does  not  predict.  Disregarding  these  peaks,  the  general  trend  for 
the  acmal  debris  flux  is  to  decrease  at  altitudes  beyond  800  to  1000  km.  The  flux  model.  Equation 
6,  does  not  predict  this  decrease  and  is  insensitive  to  changes  at  altimdes  >1000  km.  The  NASA 
90  model  appUcations  should  be  restricted  to  altitudes  <1000  km.  This  flux  model  does  not 
explicitly  account  for  variations  in  the  ballistic  coefficient  of  fragments  which  affects  their  orbital 
lifetime.  This  is  done  implicitly  in  the  solar  activity  term. 

In  rank  order,  the  relative  NASA  90  model  sensitivities  (and  comments  pertaining  thereto)  with 
respect  to  the  variables  d,  h,  i,  t,  S,  p,  and  q  are  as  follows: 

fammslSE  Comments 

•  diameter,  d  Reasonable;  however,  model  flux  becomes 

undefined  (infinite)  as  d  0.  Recommend 
limiting  model  application  to  specific  range  of 
debris  sizes,  for  example  10-^  ^  d  ^  103  cm. 

Reasonable.  Sensitivity  increases  with  growth 
rate. 

Reasonable.  Model  sensitivity  reflects  variation 
with  respect  to  debris  concentration,  4^  (i). 

Model  is  sensitive  to  choice  of  growth  (linear  or 
compounded). 

Reasonable.  Model  sensitivities  reflect 
response  or  changes  with  respect  to  t  and  S. 


•  small-particle  growth  rate,  q 

•  orbital  inclination,  i 
large-particle  growth  rate,  p 


•  time,  t 
solar  activity,  S 
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orbital  altitude,  h 


Unreasonable.  Model  predicts  essentially 
constant  flux  beyond  1000  km  with  increasing 
insensitivity  as  altitude  increases.  Limit  model 
application  to  altitudes  <1000  km. 


8.3  UNCERTAINTY  ANAT.YSTS 

Except  for  altitudes  of  500  km  (±200  km)  and  debris  sizes  >10  cm,  the  flux  model  predictions  and 
resultant  uncertainties  have  been  shown  to  be  plus  or  minus  one  order  of  magnitude  at  a  90-percent 
confidence.  The  best  measured  data  are  at  500  km  and  for  sizes  >10  cm.  The  flux  model 
predictions  at  this  altitude  and  size  range  are  within  a  factor  of  2  to  4.  The  major  contributor  to  the 
overall  model  uncertainty  is  the  uncertainty  in  the  measured  flux  environment  Reducing  only  the 
measurement  uncertainties,  Efm  and  efh,  by  an  order  of  magnitude  for  all  altitudes  and  diameters 
can  reduce  the  overall  flux  prediction  uncertainty  to  a  factor  of  2  to  4  consistent  with  500  km  and  d 
^  10  cm.  This  assumes  the  model  itself  can  be  adjusted  to  account  for  any  improved  data  or 
improved  knowledge  of  the  actual  environment  The  current  uncertainty  bound  of  at  least  plus  or 
minus  one  order  of  magnitude  (except  at  500  km  and  d  >  10  cm)  is  significant  It  is  recommended 
the  model  be  used  only  for  other  than  conceptual  design  studies. 

8.4  RESULTS/RECOMMENDATTON.S 

As  a  result  of  this  study,  the  following  conclusions  and  recommendations  are  made: 

•  The  NASA  90  model  is  an  engineering  model  for  the  prediction  of  long-term,  orbital- 
debris  environments  for  spacecraft  in  LEO  only.  It  is  an  empirical  relationship  derived 
from  curve  fits  to  the  data.  It  meets  NASA's  objective  for  a  design  model  that  is  easy 
to  use.  The  model  has  been  instrumental  in  leading  to  increased  understanding  of  the 
orbital-debris  environment 

•  It  is  recommended  this  model  only  be  used  to  establish  the  order  of  magnitude  of  the 
orbital-debris  hazard  environment  for  spacecraft  in  LEO,  specifically  h  ^  1000  km,  for 
studies  only  through  the  year  2010  and  for  present  space  operations.  The  model  should 
not  be  used  in  space  architecture  studies  where  the  space  operations  change 
significantly.  The  model  should  be  restricted  to  the  data  base,  10-4  cm  <  d  <  10^  cm. 
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•  The  collision-velocity  impact  distribution  and  direction  of  impact  algorithm  are  not 
recommended  for  quantitative  use.  The  velocity  distribution  used  in  this  model 
represents  the  expected  distribution  only  at  one  altitude  (500  km)  and  for  a  given 
snapshot  in  time.  The  estimated  relationship  for  particle  mass  rests  on  very  little  data 
and  represents  a  gross  approximation  of  the  environment. 

•  The  model  does  not  operate  with  high  accuracy  except  at  500-km  altitude  or  for  sizes 
where  d  >  10  cm.  For  all  other  regions,  the  model  uncertainty  is  at  least  one  order  of 
magnitude.  The  model  may  not  benefit  from  new  debris  measurements  to  improve  the 
model  accuracy  if  those  new  measurements  significantly  alter  the  predicted 
environment.  Presently,  the  model  does  not  reproduce  the  measured  environment 
beyond  1000  km. 

•  It  is  recommended  that  the  AF  examine  and  assess  more  sophisticated  physics-based 
computational  models  to  meet  the  goals  of  developing  an  AF,  long-term  debris  model 
that  can 

-  operate  with  high  accuracy  at  the  relevant  altitudes  (LEO)  and  orbital  parameters 
benefit  from  new  debris  measurements 

accommodate  current  and  future  AF  space  scenarios  and  architectures 

•  In  the  effort  to  develop  more  accurate  long-term,  orbital-debris  models,  ORION 
recommends  that  the  AF  pursue  a  vigorous  program  to  measure  the  debris  environment 
with  high  accuracy  in  LEO  for  all  particle  sizes. 

As  a  final  note,  the  LEO  debris  environment  is  largely  a  man-made  environment  As  obvious  as 
this  statement  may  seem,  the  development  of  any  model  which  attempts  to  predict  the  future 
consequences  of  man's  activities  must  inherently  be  uncertain. 
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